THE PERIODIC DEFOCUSING ABLOWITZ-LADIK EQUATION 
AND THE GEOMETRY OF FLOQUET CMV MATRICES 



Luen-Chau Li and Irina Nenciu 

Abstract. In this work, we show that the periodic defocusing Ablowitz-Ladik equa- 
tion can be expressed as an isospectral deformation of Floquet CMV matrices. We 
then introduce a Poisson Lie group whose underlying group is a loop group and 
show that the set of Floquet CMV matrices is a Coxeter dressing orbit of this Pois- 
son Lie group. By using the group-theoretic framework, we establish the Liouville 
integrability of the equation by constructing action-angle variables, we also solve 
the Hamiltonian equations generated by the commuting flows via Riemann-Hilbert 
factorization problems. 



1. Introduction. 

The defocusing Ablowitz-Ladik (AL) equation (a.k.a. defocusing discrete non- 
linear Schrodinger equation) is the system of differential-difference equations given 
by 

-iUn = an + l - 2an + On-l " |anP(«n-l + an+l) (1-1) 

where is a sequence of numbers inside the unit disk D. It was introduced by 
Ablowitz and Ladik in [AL] as a spatial discretization of the defocusing nonlinear 
Schrodinger equation and since then has been the subject of numerous studies. In 
particular, much attention has been focused on the inverse scattering method for 
solving the equation in the two-sided case where the index n ranges over the set 
of integers. (See [APT] and the references therein.) By contrast, the literature 
on the periodic problem is relatively sparse. (See, for example, [K], [MEKL], [N], 
[GHMT].) 

In recent years, one of the interesting developments in the arena of the defocusing 
AL equation has been the connection with the theory of orthogonal polynomials 
on the circle (OPUC) and the so-called CMV matrices [S2], [N], and our work 
here is a continuation of this development. Since we are dealing with the periodic 
defocusing AL equation here, let us begin with a set of Verblunsky coefficients 
{oj}fLo satisfying the periodicity condition aj+p = = 0,1,2,- ■■ .Without loss 
of generality, we may assume p is even. In [S2], Simon introduced the discriminant 
A(z) associated with {aj}j°,o together with the second author, they obtained 
the following involution theorem [S2], [N] 

{Aiz),Aiw)}AL = 0, {P, A(z)}al = (1.2) 
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by calculating with Wall polynomials. Here {-, -I^l is the Ablowitz-Ladik Poisson 
bracket [KM] , [S2] and P = n^=o PJ ' "^^here pj = ^\ — This prompted the 

search for an integrabic system which is related to OPUC in the same way the 
Toda lattice is related to orthogonal polynomials on the line. As it turned out, the 
sought-after integrable system is the periodic defocusing AL equation [N]. In [N], 
the Lax equations for the commuting flows were expressed in terms of the extended 
CMV matrix E with periodic Verblunsky coefficients. However, as A(z) is related 
to the Floquet CMV matrix EiK) (which is a unitary loop with spectral parameter 
K) through the characteristic polynomial det(z/ — £{h)), it is natural to ask if the 
same equations can be rewritten as isospectral deformations of £{h). As the reader 
will see in Section 2, this is indeed the case and the result is the point of departure in 
this work. More precisely, the result not only suggests that the set oipxp Floquet 
CMV matrices should have some Poisson geometric meaning, but also points to 
the linearization of such flows on geometric objects related to the Jacobi varieties 
of the underlying spectral curves. Thus our goal in this work is two-fold. First 
of all, we will link the Floquet CMV matrices to Poisson Lie groups, analogous 
to what we did in our earlier work on finite CMV matrices. (See [LI] and [KN].) 
Secondly, by using the group-theoretic framework, we will study the defocusing AL 
equation with regard to action-angle variables. We will also solve the commuting 
Hamiltonian flows via Rieniann-Hilbert factorization problems. At this juncture, let 
us mention some earlier works related to the integration of the periodic defocusing 
AL equation which is part of our second goal here. To start with, it has been known 
for quite some time that the defocusing AL equation (1.1) can be represented as 
a Lax system on a lattice (or discrete zero curvature representation), where the 
Lax operator Lj{z) associated to site j of the lattice is given by (see, for example, 
[ALl], [AL2] and [FT]) 

^^(-) = («, ?^)- M 

Therefore, in the periodic case with period p, the monodromy matrix M(z) = 
Lp-i{z) ■ ■ ■ Lo{z) undergoes an isospectral deformation which means that an equa- 
tion in Lax pair form (and different from the one we are using here) is known for the 
periodic defocusing AL equation. In [MEKL], a transformation of a natural gener- 
alization of (1.3) was discovered and the result was applied in the construction of 
finite genus solutions of a more general version of the AL equation. In particular, 
the authors in [MEKL] were able to write down the solution of the initial value 
problem for the periodic defocusing AL equation itself. On the other hand, from 
a different direction, the authors in [GHMT] considered a more general version of 
the AL hierarchy, and discussed the problem of solving the r-th AL flow when the 
initial data is the stationary solution of the p-th equation of the hierarchy. As the 
reader will see in Section 6 below, our approach in solving the commuting Hamil- 
tonian flows associated with the periodic defocusing AL equation is quite different 
from those in these earlier works. 

The paper is organized as follows. In Section 2, we begin by recalling the notion 
of CMV matrices, extended CMV matrices and Floquet CMV matrices. Then we 
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show how to rewrite the Lax equations for £ of the commuting flows associated 
with the periodic defocusing AL equation as isospectral deformations of £{h). To 
prepare for what we need in subsequent sections, we also discuss the structure of 
the powers of £{h). In Section 3, we have two main goals. The first goal is to show 
that the set of p x p Floquet CMV matrices is a symplectic leaf of a Poisson Lie 
group whose underlying group is a loop group G^. Indeed, as one would expect 
from results in [LI] and [KN] concerning the finite dimensional case, the Poisson 
structure here is also a Sklyanin structure. In fact, it is the Sklyanin structure 
associated with the Iwasawa decomposition of the loop group G^: = K^B^. 
(The decomposition was established in [GW].) However, in order to write down this 
Sklyanin structure {•, -jjit, we find it necessary to restrict ourselves to a subclass 
of functions J^(G'^) of C°°(G*). Fortunately, J'(G'^) forms an algebra of functions 
under ordinary multiplication and is closed under {-j-Ijh- Hence {-j-jjii defines a 
Poisson bracket on J^(G*). Now note that although we are dealing with a restricted 
class of functions here, it can be checked that the notion of Poisson Lie groups can 
be extended to this infinite dimensional context in a rigorous way. Moreover, we can 
check by hand that the symplectic leaves of (G* , {•, -Ijti ) are still given by the orbits 
of the dressing action. With this preparation, the technique in [LI] can be naturally 
extended to show that the set oi p x p Floquet CMV matrices is a dressing orbit 
through a Coxeter element Xf of the affine Weyl group Waff - Indeed, the induced 
Poisson structure on Kyj is a loop group analog of the Bruhat Poisson structure 
in [LW] and [Soi] and wc can show that the set oi p x p Floquet CMV matrices 
is a product of two dimensional orbits. In the rest of the section, our goal is to 
clarify the relation between the AL bracket and the Sklyanin bracket {•, -jjn , and 
to describe the Hamiltonian equations generated by the central functions on G-^j, 
thus connecting the group-theoretic framework with the equations in Section 2. 

In Section 4, wc study the analytical properties of the Bloch solution of £u = zu, 
which play an important role in subsequent sections. Since £ defines a (pcntadiag- 
onal) periodic difference operator, the seminal work of van Moerbeke and Mumford 
[MM] comes to mind. However, we note that neither £ nor its factors C and M 
in the theta- factorization of £ satisfy the genericity assumption in [MM]. So the 
analysis in this section is more delicate than the standard case [MM], [AM]. In 
Section 5, we start with a simple proof of the involution theorem in (1.2), which is 
possible because of the group-theoretic setup in Section 3. Then we proceed to con- 
struct the angle variables. To compute the Poisson brackets between the conserved 
quantities and the various quantities related to the putative angles, we make use 
of a device introduced in [DLT]. We would like to point out that in general, such 
computations could be difficult because they may require detailed information on 
the asymptotics of the normalized eigenvectors in neighborhoods of the points at 
infinity of the Riemann surface. In our case, asymptotics beyond the leading order 
are difficult to get because we are in a non-generic situation, but fortunately we 
are saved by some special structure. Finally, in Section 6 we solve the commut- 
ing Hamiltonian flows via Riemann-Hilbert factorization problems, which again are 
suggested by the group-theoretic framework. We remark that it is in this very last 
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section that we find it advantageous to think of our flows on £{h) as flows on the 

factors g'^ and g°{h) in the theta-factorization of £{h). This is precisely the reason 
why we introduce Lax systems on a period 2 lattice in Section 3. 

2. Preliminaries. 



In this section, for the convenience of the reader, we begin with some background 
material on CMV matrices and the involution theorem of Nenciu-Simon. (Good 
references are [S2] and [S3].) Then we will show how to rewrite the Lax equation in 
[N] for the periodic defocusing AL equation (in terms of the extended CMV matrix 
£) as an isospectral deformation of the Floquet CMV matrix £(h). We will also 
present a result on the structure of the powers of £{h) which we will use in Sections 
5 and 6. 

The CMV matrices are the unitary analogs of Jacobi matrices [S3] and made their 
debut in the numerical linear algebra literature. (See [B-GE] and in particular [W].) 
Subsequently, they were rediscovered by Cantero, Moral and Valazquez [CMV] in 

the context of the theory of orthogonal polynomials on the circle (OPUC). To intro- 
duce these objects, let D = {2; G C | \z\ < 1}, and let d/i be a nontrivial probability 
measure on dH, then one can produce an orthonormal basis of L'^{dn,dij) by ap- 
plying the Gram-Schmidt process to 1, -2, • ■ • . As it turns out [CMV], 
the matrix representation of the operator f{z) ^ zf{z) in L^{dO,dfj,) with respect 
to this orthonormal basis is the infinite CMV matrix 
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where aj G 



0,1, 



are the so-called Verblunsky coefficients, pj 
, and where £ = diag(^05 02,- ' -M. = diag(l, ^i, ^3, ■ 



), with 
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0,1,. 



(2.2) 



The factorization on the right hand side of (2.1) is called the t/iefa-factorization 
and lends itself to generalization. Indeed, if we now have a two-sided sequence 
{aj}'jL_^ with aj G B for all j, then we can define the extended (two-sided) 

CMV matrix £ by extending C and M. to doubly-infinite matrices in the obvious 
way. (Of course, the 1x1 block will not appear in this extension.) In this work, 
we are mainly interested in the case where the sequence {aj^JL^ of Verblunsky 
coefficients is periodic of period p and in this context, it is convenient to extend 
the one-sided sequence to a two-sided seqence {aj}^_^ satisfying the periodicity 
condition a^+p = oij for all j G Z. Thus correspondingly, we have an extended 
CMV matrix with periodic Verblunksky coefficients and such matrices have been 
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used to formulate the Lax equation for the periodic defocusing AL equation in [N] . 
Now suppose £ is an extended CMV matrix with periodic Verblunsky coefficients 
with period p. Without loss of generality, we will assume from now onwards that 
p is even (otherwise, we just replace p by 2p). Note that if T is the operator on 
Z°°(Z) defined by {Tu)j = Uj+p, then ST = T£. Therefore, if for h G dB, we define 



X^h) = {ue /°°(Z) I Tu = h-\}, 



(2.3) 



then the finite dimensional space Xf^^) is invariant under £. A basis of X(^) is given 
by the vectors 



4 = X] ^ ^^jp+k^ k = 0,---,p-l 



(2.4) 



j = -oo 



where Cj is the vector in with j-th component equal to 1 and zeros elsewhere. 

By definition, the Floquet CMV matrix £{h) is the matrix of £ \ with respect 
to the ordered basis {Sq, • • • , (5p_i), i,e. 



p-i 



£5k = ^{£{h))jk5j, A; = 0, ••• 



(2.5) 



Fom this, it is clear that the matrix of £^ \ X(^h) with respect to the same ordered 
basis is £{h)"'. Thus the entries of £{h)"' are related to those of £'^ by the formula 



{£{hr)jk = ^h-\£nj,k+ip 



(2.6) 



for < j,k < p — I. Finally, the Floquet CMV matrix £{h) also has a theta- 
factorization £{h) = g^g°{h), where 



/ = diag(6lo,6l2,--- ,9p-2), 



(2.7) 



and 







\pp-xh-^ 



7p-3 





Pp-\h\ 




(2. 



Or, 



-J 



This is of course a consequence of the factorization for the corresponding extended 
CMV matrix £. In case we want to emphasize the dependence of g*^ and g° on 
a = (q!o, • • • , Q!p-i) € ID)^, we also write 5^ = g'^ia) and g° = g°{oL)- 

Another very important notion associated with periodic Verblunsky coefficients 
is that of the discriminant introduced in [S2]: 



A(z) 



-p/2 



tr(rp(z)), zGC\{0} 



(2.9) 
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where 

is the transfer matrix. In [S2], by seeking a Poisson bracket on so that the 
modulus P = n^=o /'i generates the Aleksandrov flow, the author arrives at the 
Ablowitz-Ladik bracket (see [KM] for more general versions of this structure) 

ff f \ o- 2 f dfi df2 dfi df2\ , . 

{A,A}.. = 2.gp,(^^^-^^j^ (2.11) 

We recall the involution theorem of Nenciu-Simon which was obtained by calculat- 
ing with Wall polynomials. 

Theorem 2.1 [N], [S2]. For allz,weC\ {0}, 

{A{z),A{w)}al=0, {P,A{z)}AL=(i. (2.12) 

Hence if P-A{z) = ^^^=-^12 ^7-^"' ' the functions P, Iq, Relj, Imlj,j = 1, • • • ,p/2 — l 
Poisson commute with each other. 

This result, when combined with the proof that the functions P, /o,Re/j,Im/j, 
j = I, - ■ ■ ,p/2 — 1 are functionally independent on an open dense subset of W [S2] , 
shows that any of the functions in the above list generates a completely integrable 
Hamiltonian system. To write down the Lax pairs, the author in [N] actually 
considered a different, but equivalent set of commuting Hamiltonians, which are 
constructed from the real and imaginary parts of 



1 

Kn = -yZi£nkk, l<n<p/2-l, 



k=0 

together with Kp/2 and P. Indeed, an easy computation shows that 

{Re{Ki),aj}AL = ip]{aj-i + aj+i), {log(P), OjIal = iotj (2.13) 

for all < J < p — 1. Hence the periodic defocusing AL equation is generated by 
the Hamiltonian Re(iCi) — 21og(P). 

To relate the KnS to the Floquet CMV matrix, and to the coefficients of P- A(z), 
first recall that [S2] 

dei{zI-E{h)) = ^^Il'^jj z'^[A{z) - (/i + Zi"^)]. (2.14) 

In view of (2.13), it is clear we must consider the structure of the powers of 6. We 
will skip the proof of the following result which can be established by induction on 
n. 
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Lemma 2.2. For n > 1, {£'^)j^k is identically zero if one of the following holds: 

(a) \j — k\>2n + 1, or 

(b) j — k = 2n, where j and k are both even, or 

(c) j — k = —2n where j and k are both odd. 

As a consequence of this result, note that for n < | — 1 and for < A: < p — 1, 
we have 

{£{hr)kk = i£''Sk)k = J2 ^"'(<?")fe,.P+fc = i^^kk (2.15) 



since \k — {pq + k)\ > p > 2n + 2 for all g 7^ 0. Hence it follows that (cf. (5.4)) 

K^ = -i tr{S{hr)^, n = !,■■■ ,p/2 - 1. (2.16) 
n J\h]=i 27rz/i 

From the relation in (2.14), we see that the integral in (2.16) for n = p/2 are also 
relevant, as this is related to Iq and P. But by an induction argument similar to 
the proof of Lemma 2.2 above, we can show that 

fe+2m-l 
j=k 

and 

k-l 

{£nk,k-2m= n if A; is odd. (2.18) 

j=k—2m 

Since P = 11^=^"^ Pj for any A; G Z, it follows that 

/o ( (£P/^)kk + h-^ ■ P if /c is even 

(S(h)P/^)kk = \ J (2.19) 
K K J Jkk ^ + . p if A: is odd ^ ' 

and therefore 

tr(f(/i)P/2) = tr(fP/2) + + h-^)P. (2.20) 
Thus it follows from (2.20) and (2.15) that 

and so the set of Hamiltonians in Theorem 2.1 is equivalent to P, Rei^j, Imi^j, 
J = !,••• ,p/2-l, 

The next result gives the Lax equations of the Hamiltonian systems generated by 
the above set of functions and is the central result of [N] . We will use the following 
notation: for an infinite two-sided matrix A, n„(.A) = \Aq + v4+, where v4+ is the 
upper triangular part of A and Aq is the diagonal part. 
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Theorem 2.3. For 1 < n < p/2, 

(a) the Hamiltonian equation generated by ReK^ can be expressed as 

S=[£,iUu{S^) + i{{U^{£^)y], (2.22) 

(b) the Hamiltonian equation generated by ImKn can be expressed as 

^ = [5,n„(5-)-((n„(5"))*]. (2.23) 

Our next goal is to compute the evolution of S{h) under (2.22) and (2.23). In 
order to do this, we have to establish the following result. 

Proposition 2.4. Let n < p/2 — 1, then the structure of £{h)'^ as a Laurent 
polynomial in h is given by 

S{hT = Ao{n) + hAi{n) + /i-U-i(n), (2.24) 

where Ai{n) is strictly upper triangular and A-i{n) is strictly lower triangular. 

Proof. In order to prove (2.24), let < j, < p — 1 be two indices. Then as in 
(2.15), we have 

(^W"),. = i:^"n^"),.+,,- (2-25) 



Now, note that for |g| > 2 and any j, A;, we have the inequality 

\j - {k+pq)\>2p-\j - k\>2p- ip-l) =p+l>2n+l 

and therefore j i,_^pq = 0. Hence formula (2.24) holds for some matrices A±i{n). 
To find the structure of the matrices A±i{n), note that for j > k, we have \j — {k — 
p)\ = P+ij — k) >p>2n + l and so it follows from Lemma 2.2 that = 0. 

Consequently, 

m'')j, = {£'')j, + h-'{S-).,^^ iorj>k. (2.26) 

A similar argument shows that 

{£{hr).^ = {S-).^ + h{S-).^_^ iorj<k. (2.27) 

These last two equations then establish our claim about the triangularity of A±i{n). 
In fact, we find that 

{Ao{n)).^ = (^") for all < j, < p - 1 , (2.28) 



V ' ],k—p 


if j < k; 





if j>k, 




if j > k; 





a j <k 



iMn))^, = \ : (2.29) 

and 

iA-i{n)) .={ ^' " ■' ' ■"' (2.30) 

^ ^ ''^^ I if J < A;. 

This completes the proof. □ 

Wc are now ready to give the result alluded to above which is the point of 
departure in this work. We will make use of the projections Ilf and Il-r introduced 
in Section 3 below. (See the second paragraph of Section 3 and (3.12).) 
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Proposition 2.5. (a) If S evolves according to (2.22), then 

8{h) = [£{h),Il~,^{t£{hr)]. (2.31) 

(h) If £ evolves according to (2.23), then 

£{h) = [£{h),Il~,^{£{hr)]. (2.32) 

Proof. Let Qn{h) be the matrix of n„(£'") [ with respect to the ordered basis 
{5q, ■ ■ ■ ,6p-i). To estabUsh (2.31) and (2.32), it suffices to prove that 

Qn{h) = \{Aoin))o + (A(n))+ + /i-U_i(n), (2.33) 

where Ao{n) and A-i{n) are the matrices in (2.24), and where (^o('^))+ and 
{Ao{n))o are respectively the upper triangular part and diagonal part of Ao{n). 
But by a direct calculation and making use of Lemma 2.2, we find 

P-i fc-i 
{£-U6k = h-' ^{£-)j.,,kSj + ^{SnjkSj. (2.34) 

j=0 j=0 

Hence from (2.28) and (2.30), we conclude that the matrix of (£^")+ | ^(h) with 
respect to the ordered basis {Sq, ■ ■ ■ ,5p-i) is given by h~^A-i{n) + Ao{n)-^-. On 
the other hand, it is easy to see that the matrix of {£"')o \ X(^h) with respect to the 
same basis is (.4o(n))o. Hence (2.33) follows. □ 

We close this section with an important remark about the periodic dcfocusing 
AL equation itself. Namely, if aj = ip'j{aj-i + Oj+i), or equivalently, £{h) = 
[£{h),Ilj {i£{h))], then /3j{t) = e~'^^^aj{t) satisfies the periodic defocusing AL 
equation and vice versa. Thus in order to solve the periodic defocusing AL equa- 
tion, it suffices to solve £{h) = [£{h), H-^ {i£{h)) ], which is generated by the Hamil- 
tonian Rei^i = — Re/p/2-i(^(^))- (Compare (2.16) and (5.4).) We will solve the 
Hamiltonian equations generated by the conserved quantities in Theorem 2.1 in 
Section 6 below. 

3. Floquet CMV matrices, dressing orbits and Hamiltonian flows. 

The first goal of this section is to show that the set of Floquet CMV matrices is 
a symplcctic leaf of a Poisson Lie group whose underlying group is a loop group. 
Indeed, by following the method of investigation in [LI], we will show that there 
exist symplectic leaves of a more elementary nature in terms of which we can 
describe the collection of Floquet CMV matrices. 

As explained in Section 2 above, we can assume that p is even. Let be 
GL{p, C) considered as a real Lie group, and let K and B be respectively the unitary 
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group U{p) and the lower triangular group of p x p matrices with positive diagonal 

entries. It is well-known that admits the Iwasawa decomposition G'^ = KB. On 
the Lie algebra level, this corresponds to g"^ = t © b (with associated projections 
Hi, lib), where g'*, f. and b are, respectively, the Lie algebras of G^, K and B. 
Later on in the section, we will also need the maximal torus T of K, consisting of 
unitary diagonal matrices. 

Let G^ = C°^(5\ G^) be the smooth loop group with the (7°° topology. G^ is a 
Prechet Lie group with the Lie algebra ^ = G°°{S^,Q^). We will use the following 
nondegenerate ad-invariant pairing on q^: 

iX,Y) = Imj>^^^J.{X{h)Y{h))^. (3.1) 

As the reader will see, this choice is critical for what we have in mind. 

Following [GW], choose a symmetric weight function w : Z — y M_|_, which is 
rapidly increasing in the sense that 

lim w{n)n~^ = oo, Vs > 0. (3.2) 

n— >oo 

Also, assume that w is of non-analytic type: 

lim i(;(n)^/" = 1. (3.3) 

n— )'00 

For X G 0^ given by X{h) = E^-oo ^i^^ ^e define 

OO 

||X|U= ^ \\X,\\w{j), (3.4) 
i=— oo 

where || • || is a norm on Also, set 

(P+X)(/i) = ^X>^ {P_X){h) = Y,Xjh\ {PoX){h) = Xo. (3.5) 

Consider the Banach Lie group 

Gl = {g&&^\\\g\U<^} (3.6) 

with Lie algebra 

0«={Xe0»|||X|U<oo}. (3.7) 

From [GW], we have the Iwasawa decomposition for the loop group G^ and its Lie 
algebra 

where 

k^ = {g^Gl\g*g = I}, B^ = {g e GI \ P.g = 0, Pog e B} (3.9) 
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are Banach Lie subgroups of and 1^, are their respective Lie algebras. Denote 

by k : — > K^, h : — > the analytic maps defined by the factorization 
g = k[g)b{g)~^ , g G G*. Also, denote by Hj and Hj; the projection maps relative 

to the splitting = iw®bw Then from standard classical r- matrix theory [STSl], 
[STS2], 

Jtt=n| -Hr (3.10) 

is a solution of the modified Yang-Baxter equation (mYBE). Hence we can equip 
fl^ with the j" -bracket 

[X,YU = 1{[J^X,Y] + [X, J«y]). (3.11) 

Wc will denote the vector space g'^ equipped with the Lie bracket [•, by (05^)jtt- 
In fact, it is easy to check from (3.11) that (gj^) jtt = iwQ^w (Lie algebra antidirect 
sum). Note that explicitly, the projection maps n^^, H^^ are given by the formulas 

Ui X = p_x + h^Xq - (p_xY 

12) 

Ui^x = p+x + TibXo + {P-xy, ^' ' 

where P± are defined in (3.5). 

In order to introduce the Poisson structure on G^, it is necessary to restrict 
ourselves to a subclass of functions in G°°(G^). We say that a function ip G G°°(G*) 
is smooth at g e G* iff there exists Dip{g) G (called the right gradient of <^ at 
g) such that 



dt 



^^^^(e*^5) = (Z)^(5),X), X^qI (3.13) 



where (•, •) is the pairing in (3.1). If 9? G G°°(G^) is smooth at g for all g G G^, 
then we say it is smooth on G^. Note that the nondegeneracy of (•, •) implies that 
the map 

i--Ql^C5lr, X^{X,-) (3.14) 

is an isomorphism onto a subspace of (gj^)* which we will call the smooth part of 
(g^)*. Thus if G G°°(G^) is smooth at g iff T*rgdip{g) is in the smooth part of 
(g^)* and we can define the left gradient of such a function at g by 

j^l^^ip{ge''') = {D'^ig),X), Xeel- (3-15) 

For each g G G*, wc will denote the collection of all smooth functions at g by 
J^g(G^) and we set J'(G^) = f^g^Q' J^g{G^). With the above considerations, it 

is easy to check that T{G^) is non-empty and forms an algebra under ordinary 
multiplication of functions. 
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Proposition 3.1. (a) For (p,tp e T{G^) and g G G^, define 

{^M.n{9) = \{JKD'v{g)).D'^{g)) - i(J«(I)^(<7)),DV(ff)). (3.16) 

Then {<^, '0}jti G -^(G^) o,f^d hence {•■,■} ji, defines a Poisson bracket on J^(G*). 
(b) The Hamiltonian equation of motion generated by Lp ^ J^(G*) is given by 
9 = 9 (ni„ {D'^{g))) - (n^^ {D^{g))) g 
= {Yl~,^{D^{g)))g-g{Il-,jD'^{g))). ^ ' ' 

Proof, (a) A straight forward calculation shows that for ip,il) ^ J-(G^), D{ip, V'ljit (g) 
exists for each g € G^, and is given by 

D{^,i;U{g) =Adg[D'ip{g),D'i,{g)U+Adg^ U=o D'^pie^^^^^^^^^^ g) 



-Adg-U=o £'V(e^''(«)^'^(«)5), 



(3.18) 



'ds 

' ds 
where 

Vig) = \Adg o J« o Adg-^ - ^ J«. (3.19) 

This shows {<^, '0}jti € F{G^). To prove the second half of (a), first note that 
T^{X{h)Y{h)) G M for X,y G 1^. From this, it follows that (X,F) = 0. Conse- 
quently, fiu, is an isotropic subalgebra of g* relative to the pairing (•,•)• On the 
other hand, \i X,Y G b^,,, we have 

<p tT{X{h)Y{h))-— = tr{XoYo) G M 
J\h\=i 27rz/i 

because Xq, Iq are lower triangular with real diagonal entries. So is also an 
isotropic subalgebra of g^. Combining these two facts, we can now conclude that 
j" is skew-symmetric relative to (•,•). Finally, the Jacobi identity and the deriva- 
tion property now follow from standard calculations in [STS2] which work without 
change in our infinite dimensional context. 

(b) This derivation of the Lax equation from the Poisson structure is standard. 

□ 

Note that although we are dealing with a restricted class of functions, the notion 
of Poisson submanifolds can be defined analogously to the standard case. Now it 
is easy to check that if G J-'(G^), then so are if org and ip o Ig for all g G G^, 
where Tg and Ig denote right and left translation by g, respectively. Hence the 
notion of Poisson Lie group can be extended to this infinite dimensional context 
and (G^,{-,-}js) is a coboundary Poisson Lie group. On the infinitesimal level, 
we will call [g^, {q^)js) the tangent Lie bialgebra of (G^, {•, •} jtt ) as the map p : 
— > L{Q^,g^) {L{Q^,g^) is the space of linear maps on q^) given by p{X) = 
^{adx o J* — j" o adx) and satisfying the relation {[Y,Z]jt,X) = [Z, p{X){Y)) is 
a 1-coboundary (and hence a 1-cocycle) with respect to the adjoint representation. 
Thus in speaking of a Lie bialgebra here, the underlying vector spaces of the pair 
of Lie algebras involved are only required to be in duality with respect to (•, •) and 
this is what we will continue to do. 
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Corollary 3.2. (a) K^, is a Poisson Lie subgroup of {G^, {•, - jjit) in the sense that 
Kw is a Lie subgroup of which is also a Poisson submanifold o/ (G* , {•, •} jtt ). 
Moreover, the tangent Lie bialgebra {iw,{9^)js/^w) '^f (where is defined 
relative to (•, ■)) is isomorphic to (t^u, b~) where b~ is but equipped with the — 
bracket. 

(b) The underlying group of the Poisson group {G^)js dual to (G^, {■, -Ijii) consists 
of equipped with the multiplication 

g*h = k{g)hb{g)-K (3.20) 



Proof, (a) To show that Ky^ is a Poisson Lie subgroup, it is enough to check that 
Ku, is a Poisson submanifold of (G^, {•, -jjii) and this can be done by using the 
expression for the Hamiltonian vector field in Proposition 3.1(b). To show that the 
tangent Lie bialgebra (^u., (0*)jt /t^) of /v,^, is isomorphic to (t^,b~), first note 
that 1^ =tu,- Hence we have X + 1^ = Hg^X + 1^ for all X e {g^)ji. Therefore, 
the induced Lie bracket on (g*) js /i^ is given by 

Conseqently, the map X + X is an isomorphism, when b^ is equipped 

with the — bracket. 

(b) The formula is a consequence of the fact that 0^ = and can be verified 
easily. □ 

Remark 3.3. In view of the second half of Corollary 3.2 (a), the induced structure 
on Kw is the loop group analog of the Bruhat Poisson structure in [LW] and [Soi]. 

We next turn to the description of the symplectic leaves of the Sklyanin structure 
in (3.16). Unfortunately, we cannot assume the general results in [STS2] and [LW] 
apply to our case without some verification, because the analysis in these works is 
for finite dimensional Poisson Lie groups. In this regard, let us also remark that 
as far as we know, the integrability of the characteristic distribution of a Poisson 
structure on an infinite dimensional manifold is by no means automatic because an 
analog of the Stefan-Sussmann result [St, Su] is not available. In the following, we 
will check things by hand. So let us define 

S, = {X^{x)\^eT,{Gl)} (3.21) 

for each x G G^. Then the characteristic distribution of the Poisson bracket {•, -jjtt 
is given by 5 = Li^^Qw S^- On the other hand, if (G^)}it denotes the identity 

component of {G^)j», the right dressing action [STS2] of {G^)j» on is defined 
by the formula 

^n(x) = k(q)~^ xk(x~^ qx) 

_i _i (3-22) 
= Kg) xb{x gx) 
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and the infinitesimal generator of this action corresponding to ^ G is the vector 
field on G^, where 

VJ 

Cg« (^) = ^xJ^x-'^x) - (3.23) 

For each x G G?,, let be the subspace of T^G^ spanned by the vectors ^^=,3 (x). 
Then F = U F^ is an integrable generalized distribution whose leaves are the 
orbits of the dressing action 

Proposition 3.4. For each x G G^, T^O = S^, where O is the orbit of the dressing 
action containing x. Hence the characteristic distribution S is integrable and the 
leaves of this distribution are given by the orbits of the dressing action 

Proof. From (3.7), 

X^{x) = ^xJ^{x-^Dip{x)x) - ^j\Dip{x))x. (3.24) 

To show that = F^, it suffices to show that for each ^ G g^, there exists 
ip G J^xid^) such that Dip{x) = ^. To do so, we use the fact that the exponential map 
exp : gj^ — > G^ is a diffeomorphism of a neighborhood ?7 of onto a neighborhood 
V of the identity element of [GW]. Clearly, Vx is a neighborhood of x and in 
this neighborhood, define (p^{g) = {^,log{gx~^)) and extend this to a function 
(f^ G C°°{Gl). Then ip^ G T^^iG^) with D(p^{x) = ^. This completes the proof. □ 

Let O be a dressing orbit, as in the proposition above. We next show that 
O is a symplcctic leaf of {G^^,, {■, ■} jt), i-c, there exists a weak (resp. strong) 
symplectic form (see, for example, [OR] for such matters) uiq on O consistent 
with the Poisson bracket {•, -Ijh in the case when O is infinite (resp. finite) di- 
mensional. For this purpose, let (gl^)^ denote the smooth part of (g*)* and let 
iT*Gl)s = U,eGi^^.-(0^)5- Also, denote by ^« : iT*Gl)s TGl the bun- 
die map corresponding to {•, -jjs and let j be the left inverse of i : g^ — > (g^)*. 
For each x G O, we define a skew-symmetric bilinear form uj^ on Sx by the formula 

cOx{Tr\x){a),J{xm) 

= {Txr,-^7r\x){a),j{T:r,P)) (3.25) 

= -(i(r;r,a),r,r,-.7r«(x)(/3)). 

Clearly, the value of the above expression depends only on the values of 7r^{x){a) 
and 7r''(x)(/3). Thus is a well-defined skew-symmetric bilinear form on Sx- Now 
suppose uJx{iT'^{x){a),7r^{x){j3)) = for all 7r''(x)(/3). Then from (3.25) and the non- 
degeneracy of (•, •), it follows that j{T*rxOt) = which in term implies 7r''(x)(a) = 0. 
So this establishes the nondegeneracy of ojx- Thus there exists a 2-form ojo on O 
such that 0Jo{x) = uix for each x € O. Now the argument that uq is differentiable 
and closed follows as in the finite dimensional case in [Ko] . Conseqently, uo defines 
a Poisson structure {•, - jo on O and we have 

{ip,i^}j,\0 = {ip\0,i;\0}o (3.26) 
for all ip,ip G J^(G^). Hence we have established the following result. 
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Proposition 3.5. The symplectic leaves of (G^, {■, -Ijs) are given by the orbits of 
the dressing action 

In our next step, we will make the eomiection between certain symplectic leaves 
of (GJ^, {•, •} j3 ) (these are also those of the Poisson Lie subgroup K^) and the 
Floquet CMV matrices. We begin with some notations. As in [LI], we will denote 
by g^ any p x p block diagonal matrices with 2x2 diagonal blocks of the form 



a p 
p —a 



with a G 



and p = ^/l — 



a\ 



(3.27) 



We will denote the collection of such matrices by T^. On the other hand, we will 
denote by g°{h) any loops in of the form 



/ a 








\ph-^ 



ph\ 






-aj 



(3.28) 



where a G D, /o = y^l — |aP and 6 is a (p — 2) x (p — 2) block diagonal matrix 
with 2x2 blocks of the same kind as in g^. We will denote the collection of such 
unitary loops by T°. Clearly, for given g^ G T*^ and g" E T°, the product g'^g°{h) 
is a Floquet CMV matrix. Indeed, the map 



m\T^ xT° :T'' xT° — y {pxp Floquet CMV matrices} 



(3.29) 



is a diffeomorphism, where m : x Ky^ — > is the multiplication map of 
the Poisson Lie subgroup K^j,. Finally, we will denote the dressing orbit through 

X G by O^. 

In analogy to formula (2.21) in [LI], we introduce the following special Floquet 
CMV matrix 

Xf{h)=x}x"f{h) 



corresponding to 
In other words, 
and 



a = (0,0,-- - ,0). 

x'f = diag(w*, w*, • • 



( 










w* 






h\ 






o/ 



(3.30) 
(3.31) 
(3.32) 

(3.33) 
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where 



w 



1 

1 



(3.34) 



These matrices are elements of the affine Weyl group Waff = Wt<T [PS], where 
W = N{T)/T is the Weyl group of K, and T is the lattice of homomorphisms 
— > T. Indeed, if f) is the Cartan subalgebra of Q consisting of diagonal matrices, 
Xi — Xj, i j are the roots corresponding to the pair (g, f]), then in terms of the 
simple roots a, = Aj — Aj+i, i = 1, ■ ■ ■ ,p — 1 and the highest root 6* = Ai — Ap, we 
have 

(3.35) 



X4 



and 



X'f{h) = Woc.,Wa^ ■ ■ ■ Wap_^ exp(/l(£'oo - Ep-i^p-i))w0. 



(3.36) 



Here for each j = 1, ■ ■ ■ ,p — 1, w^j = diag{Ij-i,w* , Ip-j-i) is the element in W 
which corresponds to the simple reflection Sq, , while 



We 



/O 





Vi 



1\ 




0/ 



(3.37) 



is the element in W which corresponds to the reflection sq. Finally, the element 
exp(/i(£'oo ~ -^p-i,p-i)) is in T, where Ejj denote the diagonal matrix with a 1 in 
the position and zeros elsewhere. But now recall that there is an additonal 

element ao = 6 — 9 in the simple system of affine roots in addition to those given 
by the extensions of the a^'s, j = 1, • • • ,p — I. (see, for example, [Mac] for details 
on affine Lie algebras). With ao, we can interpret the product w^^ = cxp(/7,(£'oo — 
Ep-i^p-i))w0 as corresponding to Sq,q. Therefore, by (3.35) and (3.36), we conclude 
that the CMV matrix Xf{h) introduced above is a Coxeter element of the affine 
Weyl group. With this background, we are now ready to accomplish our first goal 
of this section. 



(a) Oxj 



X^f 



n x'} 



Theorem 3.6. 

(h) o^o = 

(c) Equip X with the product structure, then Ox'j x Ox° is a symplectic leaf 

of X K^u. Moreover, the collection of p x p Floquet CMV matrices is the image 
of Oxj X Ox° under the Poisson automorphism m \ Ox<=^ x Ox° ■ C^j x Ox° — > 

{p X p Floquet CMV matrices }, where m is the multiplication map of Ky,. Hence 
{pxp Floquet CMV matrices } = = Ky, fl B^^ x / By, , a Coxeter dressing orbit. 

Proof, (a) Take an arbitrary element 



a = k{g)-^x}k{{x})-'gx}] 
= b{g)-'x}b{{x})-'gx}) 



(3.38) 
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in the dressing orbit through x^. Prom the first hne of the above expression, it 
is clear that a is unitary. On the other hand, it follows from the second line of 
the above expression that P_a = and that Poa is block lower triangular with 
2x2 blocks on the diagonal. Moreover, from the fact that the elements in B have 
positive diagonal entries, it follows that each of the 2x2 blocks on the main diagonal 
of PqU has the following properties: (i) the entry in the upper right hand corner 
is positive, (ii) the determinant is negative (since det{w*) = —1). But P_a = 
implies P-a~^ = 0. As P-a~^ = P-a* = (P+a)*, we conclude that P+a = 
and hence a = Pga. But then [a*)~^ = ((Poa)*)""'' is upper block triangular with 
diagonal blocks having the same properties. Since a = (a*)~^, it follows that a 
must be block diagonal, i.e.. 



a = diag((?!)o, (t>2,--- , 4>p-2) (3.39) 

where for each j, (f)2j is a unitary matrix with a positive entry in the upper right 
hand corner and whose determinant is —1. Consequently, (f)2j must be of the form 

' \p2j -a2jj 

for some a2j G B, where p2j = (1 — |Q!2j|^)^- Hence we have shown that O^j C T^. 
The reverse inclusion T'^ C Oxj follows exactly as in the proof of the corresponding 
assertion in Theorem 2.4 (a) of [LI], 
(b) Take an arbitrary element 

b = k{9)-'x"fk{{x})-'9x°f) 
= bi9)-^x'}b{ix'})-^gx}) 

in the dressing orbit through Xj. From the first line of the above expression, b is 
unitary. On the other hand, it follows from the second line of the same expression 
that 

oo 

b{h) = J2 hh'. (3.42) 

i=-i 

Now, 

b-' = {b{{x})-'gxl))-\x"^rb{g) (3.43) 
since x°^ is unitary. Prom this, it follows that 

{P+{h-^h)Y = P-{(hb)*) = P-{hb-^) = 0. (3.44) 



Therefore, when we combine this with (3.42), we conclude that 



b{h) = b_ih~^ +bo + bih. 



(3.45) 
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Considering the coefficient oi h ^ in the second hne of (3.41), we see that 

6_i = (6(5)-^)o (6((a;p- VP)o > 

as the first and last factors in the second hne of (3.39) only contain nonnegative 
powers of h. Since {x°i)-i has only one nonzero entry in its bottom left corner, and 
it is multiplied by lower triangular matrices with positive diagonal entries, we see 
that all the entries of 6_i are zero, except (6_i)p_i_o = 7 > 0. 

Further note that, if x{h) is unitary for all h E S^, then an element y{h) G if 
and only if y{h)* G Ox*- Indeed, 



y{h) = k{g)-'xk{x-'gx) ^ y{hy = k{g')-^x*k{xg' x*), 

where g' = x~^gx. In our case, (xj)* = x°j, and hence 

b{h)* =blh-^ + b*o + b*_ih eOxo. 

So we see, by the above argument applied to bl , that bi must have only one nonzero 
entry, in its upper right-hand corner: (6i)o,p-i = 6 > 0. 

It remains to understand bo. Recall that b{h) is unitary and expand the right- 
hand side of 1 = b{h)* ■ b{h) = b{h) ■ b{h)* in powers of h. Given the structure of 
the matrices b±i, it is immediate that the coefficients of /i^^ are both zero. The 
coefficients of must also be zero. This translates into 



and 



6^50 + blb-i = 



6_i6* + bobl = 0. 



Taking into account the shape of and writing down explicitly these relations 
leads to the fact that 



where 



0-th row of 6o = (r/,0... ,0), 
[p —1) — th row of 6o = (0, . . . , 0, i^), 

Sr] + ^1? = 0. 



(3.46) 
(3.47) 

(3.48) 



Therefore, when we combine our analysis above on together with (3.46) and 
(3.47), our conclusion is that b{h) looks like 











\-ih-^ 



5h\ 






V J 



(3.49) 
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where 6 is a (p — 2) x (p — 2) matrix which is independent of h. But a unitary 

matrix can only be a direct sum of unitary matrices, and b{h) is unitary. Unitarity 
of the 2x2 piece at /i = 1, together with the fact that 7, (5 > 0, means that 



7/1"^ u J \Pp-ih 



1 „ i-l . ..-1 . ) (3-50) 



for some complex number ap-i G B and Pp-i = \/^ — |Q!p_i p. 

The last step requires that we understand the shape of h. This is achieved by 
writing the coefficient of /i° in 6(/i) and the second line of (3.41). Again using 
the fact that the coefficients of /i° in the external factors on the right-hand side 
of (3.41) are lower triangular with positive diagonal entries, and keeping track of 
all the (potentially) nonzero entries, we obtain that 6 is a unitary matrix which 
is block-lower triangular, each 2x2 diagonal block having a positive entry in its 
upper right-hand corner. Hence the same argument as in part (a) of this theorem 
shows that 

6 = diag(^i,03,...), (3.51) 



(3.52) 



where 

' "fe Pk 
Pk -Oik 

Inserting this into (3.49) and using (3.50) shows that h G T° . 

The proof is finished once we prove that T° C O^-". This is done as in the 
previous cases. More precisely, if 5 = 6*1 © ^3 © • • • ® 9p-i{h) is an element of T°, 
choose a matrix g = li® ■ ■ ■ ® lp-3 © /p_i(/i), where, for < j < {p — 4)/2, 

is a 2 x 2 block situated between the rows and columns 2j + 1 and 2j -|- 2, and 



Pp. 



1 



Note that hj+i are all independent of h and lower triangular with positive diagonal 
entries. Further note that 

In particular, this implies that g G and hence 

k{g)-^x'}k{ix})-'gx'}) = k{gx"f). (3.56) 
Furthermore, we have the factorizations 



(3.57) 
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for < j < (p-4)/2, and 

In other words 9 ■ x'j = b ■ g, where g is in B^j for the same reason as g is. We 
therefore conclude that 

kig)-'x'}ki{x;)-'gx'}) = kigx}) = b (3.59) 

is indeed an element of O-r" ■ 

(c) It is easy to see that O^. = O^" ■ O^" ■ The rest of the assertion is clear from 
what we have already done. □ 

Remark 3.7. By essentially following the same argument, we can also establish the 
following fact: 

{pxp Floquet CMV matrices } = O^^^ ■ ■ ■ 0„^^_^ -O^^^ ■ ■ ■ 0^^^_^ -O^^^ , (3.60) 

where each orbit on the right hand side is two-dimensional. 

In the next result, we clarify the relation between the Ablowitz-Ladik bracket in 
(2.11) and the Sklyanin bracket {-j-jju. 

Theorem 3.8. The map 

BP — > G^, 

(3.61) 

a = {ao,--- ,ap-i) ^ £{h) =/(a)5°(a)(M 

is a Poisson embedding, when is equipped with the Ablowitz-Ladik bracket, and 
is equipped with Sklyanin structure {•,-}jtt. 

Proof. As the multiplication map of is a Poisson map, it is enough to show 
that the map a ^ (^^(a), ff°(a)) is Poisson, when x G^ is equipped with the 
product structure. For this purpose, denote by Eab the pxp matrix whose (a, b) 
entry is equal to 1 and whose other entries are zero. For j even, j € {0, • • • ,p — 1} 
and for I odd, I € {O, - ■ ■ ,p — 1}, introduce the following functions on G^ x G^: 



Fj(A,B) = Im £^_tr (A(h)E,j) Oj{A,B) = Re jf tr (A{h)E,,) 

F,(A. B) = Im £ iB(h)Eu) ^. GM, B) = Re £^ {B(k)E„) ^. 

(3.62) 

Then we have 

Fj {g^ ,g°) = -Im aj , Gj {g^ ,g°) = Re aj 
Fi{g',g°) = -Imai, Gi{g',g'') = Reai. 
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In view of this, it suffices to compute the Poisson brackets of these functions at 
(5*^) 5°) £ X G^. For a function ip on x which is smooth in both variables, 
we denote by Di(p{A,B) (resp. D^(p{A,B)) its right gradient (resp. left gradient) 
with respect to the z-th variable, i = 1,2. Then we have 

D[F,{g\g°) = E,,g^ = a,E,, + p,E^j+^, (3.64) 
D2Fj{g%g°) = 0, D'^F^ig^ g") = 0, 

and similarly, 

D2Fi{g\g°) = g°{h)Eu = aiEu + piEi+i^i, 

D!,Fi{g^,g") = Eug"(h) = aiEu + piEij+i, I odd, l<l<p-3, 

D'2Fp_i{g'',g") = g°{h)Ep_i^p_i = pp_ihEo^p_i + ap-iEp_i^p_i, (3.65) 

D2Fp-i{g^ , g°) = Ep_i^p_ig"{h) = Pp-ih ^£^p_i,o + 

D^Fi{g-,g'')=Q, D[Fi{g^ , g") = Q, / odd, 1 < / < p - 1. 

On the other hand, it is clear that 

AG',•(<7^<7°) = iD,Fjig-,g''),D[G,{g-,g") = iD[Fj{g'^ , g°), 
D2G,(/,5°) = 0, L>^G,(/,5°) = 0, 

and similarly, 

D2Gi{g\g") = iD2Fi{g^,g°),D'2Gi{g%g") = iD'2Fi{g',g°) 
D^Gi{g-,g°) = 0, D[Gi{g'^ ,g°) = 



(3.66) 



(3.67) 



In what follows, the indices j and k arc even, while the indices I and m are odd, 
and all indices are from {0, • • • ,p — 1}. Let {•, denote the product structure on 
G^ X G^, then it is immediate from the definition of {•,•}* that 

{Fj,Fi}, = 0, {Gj,Gi}, = 0, {Fj,Gi}, = 0, {Gj,Fi}, = 0. (3.68) 

Now, from (3.64), DiFj{g^,g°) and D[Fj{g^,g") are constant loops, so it follows 
from (3.12) and equation (2.6) of [LI] that 

J«(L>iF,(/,5'')) = -ajEjj - 

j\D[Fj{g^,g")) = p^Ej,j+^ - 2pjEj+,,j - ajEjj. 

Therefore, on using (3.69) and (3.64), we find for j ^ k that 

tr {J«{D[F,{g\g"))D[F,{g',g'')) 
={D[Fk{g^g"))kk{JHD[Fj{g',g")%k + {D[Fk{g^g")%^^^ 
= 0, 

(3.70) 
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and similarly, 

tr {.P{DrF,{g',g°))DrFk{g',g°)) 

HDiFk{g',g"))kk{JKD,Fj{g\g'^))kk + {DiFk{g\g°))k+^^^^ 
= 0. 

(3.71) 

Thus it follows from (3.70), (3.71) and (3.66) that 

{F,-,Fj,(/,5°) = 0, {G,,Gfc},(/,3°) = 0, {F„Gk}.ig',g°) = 0, for j / fc. 

(3.72) 

In a similar fashion, by using (3.64), (3.66). and (3.69), we find 

{F„G,}.(ff^ff°) 

= hr {zj\D[F,ig^,g''))D[F,{g^,g")) - ^tr {i J\D^F,{g^ , g'^))D,F,ig^ , g'')) 

= -pI 

(3.73) 

Analogously, for the odd indices Z, m, with Z, m G {0, • • • ,p — 3}, we have 

{Fi,F^}.(/,3«) = 0, {G/,G^}.(/,5°) = 0, {F;,G^}.(/,5°) = 0, for I / m. 

(3.74) 

Also, 

{Fi,GJ*(5%5'') = -pf- (3.75) 

Next, we consider brackets with the quantities Fp-i and Gp-i. To this end, we note 
the formulas 

J^{D2Fp_i{g^,g")) = -ap_iEp_i^p_i - Pp-ihEo^p_i, 

-1 (o.7d) 
J''(£>2-fl>-l(5^5°)) = Pp-ih -Ep_i,o - ap_iEp_i^p_i - 2pp_ihEo,p-i. 

Therefore, if I is odd, Z < p — 3, a calculation similar to (3.70) and (3.71) above 
shows that 

{Fp_i,Fj,(/,5°) = 0, {Gp_l,GJ.(5^5°) = 0, {Fp_^,GiUg-,g°) = 0. (3-77) 
Now, by (3.76), (3.67) and (3.65), we have 

tr (jt(D2Fp-i(/,5°))I?2Gp_i(/,5°)) = (3.78) 

while 

tr (J«(Z)^Fp_i(/,5°))Z)^Gp_i(/,<7°)) = -2ipl_, - (3.79) 
Consequently, 

{Fp_i,Gp_i},(/,5°) = -pl-r- (3.80) 
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Assembling the calculations, we conclude that for all a, 6 G {0, • • • ,p— 1}, we have 
{Fa,F,},{g^,g'^) = 0, {Ga,Gb}*{g' , g") = 0, 5°) = 0, if a 7^ 6, 

{Fa,GaUf,9") = -pl 

(3.81) 

So finally, wc obtain the following Poisson bracket relations 
{Ga - iGa,Gb - iFkUg^g") = 0, {G„ + + 5°) = 0, if a / 6, 

{Ga - iFa,Gb + iFbW,g°) = 2i6abPl, 

(3.82) 

as desired. 

□ 

Finally we describe the Hamiltonian equations generated by central functions on 
G^ in the above framework. We also introduce the kind of equations which we will 
need to use in Section 6 below. 

Proposition 3.9. (a) The Hamiltonian equation of motion generated by a central 
function ip on G^ is given by the Lax equation 

9 = g (n^^ (D^ig))) - (n^^ {D^{g)))g 
= iU~,jD^{g)))g-giU-,jD^{g))) ^ ' > 

(b) Consider (G^, {•, - jjit) and equip the group G^ x G^ with the product Poisson 
structure. If is a central function on G^., then the Lax system 

91 =91 {U~^^{D^{g2gi))) - {U'^jDip{gig2)))gi 
= (n^^ iDip{grg2))) 9i - 9i (n^,^, {D^{g29i))) , 

92 = 92 {U~^^{Dcp{gig2))) - {U~^jDip{g2gi))) 52 

= {U~,jD^{g2gi))) 92 -92 {li~^jD^{9ig2))) • 

is the Hamiltonian equation on G^ x generated, by H^p{gi,g2) = Lp{gig2)- More- 
over, under the Hamiltonian flow defined by (3.84), 9 = 9i92 evolves according to 
(3.83). 

(c) Ifki{t), bi(t), i = 1,2 are the solutions of the factorization problems 

gtDcp(si (0)32(0)) ^ ki{t)b^\t), e'D^^a^{o)g,m ^ ^2(t)62 '(t), (3.85) 

where ki (t) G , 6j (t) G , then the flow defined by (3.84) given by 
9iit) = k^{t)-'g,iO)k2it) = b,it)-' g^iQ)b2it) , 
92{t) = k2it)-'g2iO)k,it) = b2it)-'g2iQ)b,{t). 

To conclude this section, we remark that equations of the type in (3.84) are a 
special case of so-called Lax systems on a periodic lattice or difference Lax equations 
and we refer the reader to [STS2] and [LP] for the general theory. In Section 6 below, 
we will show how to solve the factorization problems for the flows generated by the 
commuting integrals of the periodic defocTising Ablowitz-Ladik equation by means 
of Riemann theta functions associated with a hyperelliptic curve. 
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4. Analytical properties of the Bloch solution. 

For any z € C, consider the equation 

Su = zu, (4-1) 

where £ is the extended CMV matrix with periodic Vcrbhmsky coefficients with 
period p, as in Section 2. Since £ admits a ^-factorization £ = CM [82], as in the 
one-sided case, it follows that (4.1) is equivalent to 

Mu = C*u. (4.2) 

In terms of the components of u and the entries of £ and M, (4.2) gives the three- 
term recurrence relations 

p2j-lU2j-l - a2j-lU2j = z{a2jU2j + p2jU2j+l) 
Oi2j+lU2j+l + P2j+lU2j+2 = z{p2jU2j - a2jU2j+l) 

for all J G Z. Due to the equivalent form in (4.2), the space of solutions of (4.1) is 

two dimensional. Indeed, it is clear from (4.3) that for given values of n_i and uq, 
we can determine all other values of by recursion. For our analysis, we will fix 
a basis with the following initial conditions: 



<P-i{z) = l, Mz) = o, 
V'-i(z) = o, Mz) = i- 

By using the first relation in (4.3) corresponding to j = 0, we have 



(4.4) 



Mz) = (4.5) 

zpo 

In general, an easy induction using (4.3) gives the following result. 
Proposition 4.1. For all j > 1, 



<^ .(^) = ^0Pp-i_^,_i a2j-iPp-i 1^ 

P0---P2j-1 P0---P2j-lZl 

aoa2jPp-i ,_i Pp-i 1 

(f)2j + l{z) = Z^ +■■■ + 

P0---P2j PQ---p2jZ^ + ^ 



(4.6) 



In particular, 



Similarly, we have 



PO • • • Pp-2 PO • • • Pp-2 Z'pI'^ 



V'i(^) = , (4.8) 

Po Po z 

and by induction, we obtain the following analog of Proposition 4.1. 
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Proposition 4.2. For all j > 1, 



i^2j{z) = + • • • + — ^ -, 

P0---P2j-1 P0---p2j-lZ^ 

I f \ ^■'ij 7 OLy-\ 1 

V'2i + l(2) = -- J ~J 



(4.9) 



In 'particular, 



^P-i(-) = ~-^^^z^l--^ -^^^-i,, (4.10) 



Po • • • Pp-2 Po • • • Pp-2 



la iP 1 

V;p(z) = + . . . + i-. (4.11) 

Now by the periodicity of the Verblunsky coefficients, we have 

(^,+p(z) Vi+p(^)) = (0,(-^) V'i(^))M(z) (4.12) 

for all J, where 



(4.13) 



is the monodromy matrix. 

Proposition 4.3. For all z, detM(2) = 1. 

Proof. Let Wj{z) = pj{(f)j{z)'tpj-^i{z) — (f)j-^-i{z)'tpj{z)). Then from the first relation 
in (4.3), we have 

W2,-iiz) = -zW2j{z) (4.14) 
for all j. Similarly, from the second relation in (4.3), we find that 

W2j+iiz) = -zW2jiz) (4.15) 

for all j. As the right hand sides of (4.14) and (4.15) are equal, it follows that 
W2j-i{z) is independent of j and consequently Wp-i{z) = W-i{z) from which the 
assertion follows. □ 

From Proposition 4.3, the eigenvalues of the monodromy matrix (i.e., the Floquet 
multipliers) are the roots of the characteristic polynomial 

-tTM{z)h + l = (}. (4.16) 

If T : Z°°(Z) — )■ /°°(Z) denote the shift operator defined by {Tu)j = Uj+p, then 
the unique solution of the problem 

£f = zf, Tf = h-'f, fp-i = l (4.17) 
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is called the Bloch solution and finding this solution is equivalent to considering 
the spectrum of the corresponding Floquet CMV matrix £{h), 

£{h)v = zv (4.18) 

where 

/ /o \ 



fp-2 

V 1 / 

Hence the ordered pair {z, h) in (4.17) must obey the equation 



(4.19) 



detizi - £{h)) = ^YIpjJ z^[Aiz)-{h + h-^)]=0, (4.20) 
where the discriminant A(2;) is related to the transfer matrix 

by the formula [S2] 

Aiz) = z-P/^tiTpiz). (4.22) 

By comparing (4.16) and (4.20), we therefore conclude that trM(2;) = A(z) and 
this relates the multiplier curve and the spectral curve. We will make the genericity 
assumption 

{GA)i the roots of P(z) = (n?=o P^)^ (A^z)^ zP - AzP) = UZii^-^^i) are distinct. 

Then it is straightforward to check that the affine curve as defined by the equation 

I{h,z) :=h -detizi -S{h)) = (4.23) 

is smooth with branch points located at Ai,--- ,X2p- We will denote by C the 
hyperelliptic Riemann surface of genus g = p — 1 corresponding to this afhne curve. 
In order to find the divisor structure of j = 0, • • • ,p — 2 on C, let 

(z) = Q+ + Q- - P+- P- (4.24) 

where P+,(5+ are on the + sheet and P-,Q- are on the — sheet of C. (The it- 
sheets correspond to the choice of sign in front of the radical in the first line of 
(4.25).) Solving for h in terms of of z from (4.20), we find 

zP/^A{z) ± ^/A{z)^zP - AzP 



h{z) 

A(z) ± zPl-^A^z) \\ - + • • • 1 

for z near oo, 



(4.25) 



+ • • • for z near oo on the + sheet, 

11^=0 Pi 



J = 

p-l 

+ ■ ■ ■ for z near oo on the — sheet. 
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On the other hand, it is clear from (5.1) that P z^l"^ l^{z) ~ 1 as 2; — )■ 0. Therefore, 
it follows from the first line of (4.25) that h{z) ^ p as z — >■ on the + sheet, 
while h{z) ~ P z^l"^ as 2; on the — sheet. Thus we have 

{K) = + |P_ - |q+ + |q_. (4.26) 

Proposition 4.4. For each < j < p — 2, 

f,{P) = h{P)<i>MP)) + ^^^^^|^^y^^^(^(^))' P ^ (4.27) 

Proof. Since cf) and if) form a basis of solutions of the equation £u = zu, we must 
have fj = ci<pj + C2i{jj for some constants ci and C2. Putting j = —1 and in the 
above expression and using the initial conditions in (4.4), we find fj = h(pj + /oV'j- 
As the vector with components h and /o is an eigenvector of the monodromy matrix 
M{z) with eigenvalue h~^, wc find that (pp-ih + V'p-i/o = 1- Solving for /o from 
this expression, we obtain the desired expression for fj . □ 

We next make the following assumption. 

{GA)2 a,-7^0for j = 0,--- 

Note that in particular, ap_2 7^ and so the degree of the polynomial z^/^^-'p-i (•?) 
isp— 1. The roots of this polynomial will be denoted hy Zk, k = 1, ■ ■ ■ ,p. (In general, 
the Zfc's are not necessarily all distinct.) 

Proposition 4.5. For each < j < p — 2, fj{P) is a single-valued meromorphic 
function on the Riemann surface C . On the finite part of C away from Q+, fj{P) 

has at worst poles at the points = (</'p-i(zj[j;), z^), where ipp-i{zk) = for k = 
1, • • • ,p—l. Moreover, the Zk 's coincide with the eigenvalues of the Dirichlet problem 

£u = zu, u-i = 0,Up-i = 0. (4.28) 

Equivalently, the z^ 's are the zeros of the equation 

det {z{fy - f{h)) = 0, (4.29) 

where [g^)* and g°{h) are (p — 1) x (p — 1) matrices obtained from {g")* and g°{h) 
by removing their last row and last column. 

Proof. For each j, it follows from (4.27) that fj{P) is meromorphic on C. Let 
us consider a point Zk where ^p-i{zk) = 0. From (4.13), we see that at such a 
point, the Floquet multipliers are given by (pp-i{zk) and (f)p-i{zk)~^ and hence 
{(f)p-i{zk), Zk) and {4>p-i{zk)~^ , Zk) are points on C. Clearly, 1 — h4>p-i{z) van- 
ishes at {4'p-i{zk)~^ , Zk). Therefore, provided that (pp-i{zk) ^ ±1 or equivalently, 
A[zk) / ±2, fj{P) has a pole at Pk = {4'p-i{zk), Zk). Since tp-i = 0, we see that 
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the solutions of = coincide with the eigenvalues of the Dirichlet problem 

in (4.28). On the other hand, observe that if il)p-i{zk) = 0, then 



diag(-Q;p_ 1,6*1, • • • ,Op--i) 



: = Zk diag(6io, • • • , 6i*_4, ap_2 

\'^p-2{Zk) ) 




from the connection with (4.28) where iI)q = 1. As the matrix on the left hand side 

of the above formula is g°{h)^ while the one on the right hand side is (5^)*, the last 
assertion in the proposition follows. □ 

Remark 4-6- The Dirichlet eigenvalues above should not be confused with the 
Dirichlet data in Chapter 11 of Simon [S2]. While the latter quantities always 
reside on the unit circle and the number of such quantities is equal to p, it is not 
the case for the z^s, as is evident from the relation rft-i Zk = —^r^- 

In what follows, we will denote by h^{z) (resp. f^{z)) the values of the function 
h{P) (resp. fj{P)) on the ± sheets of the Riemann surface. 

Proposition 4.7. For j = 0, ■ ■ ■ ,p/2 — 1, we have 



/2i(^) ~ _ ^ Z- 



^ -(p/2-i-l) 
i=2j+l 



(4.30) 



as z 00. Hence j2j{P) o-nd /2j+i(P) have zeros of order p/2 — j — 1 at P_. 

Proof. Consider first the even case. By using (4.25) and Propositions 4.1 and 4.2, 
we have 

as z ^ 00. Similary, 



ft) 



h-{z)M^) ' (4-32) 



as 2; ^ 00. Therefore, on comparing (4.31) and (4.32), the assertion for the even 
case follows. We will skip the details for the odd case as it proceeds in the same 
way. 

□ 
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To investigate the behaviour of fj~{z) as z — )■ cxd, we will first establish an identity 
for the product fj~{z)f^{z). To this end, observe that 



1 - h^(j)p-i{z) 
1 



(2 - cf)p_i{z)i4>p-iiz)+i;p{z))) + </.p_i(z)VA(z)2-4 



Therefore, by a direct multiplication and using (4.31), we find 

j+,p.j-,p. _ -tljj{(j>j(l)p-i + ipj(j)p) + (pjjtpjipp + (pjtpp-i) 
^ ^ tpp-i 



(4.33) 



(4.34) 



where on the right hand side, we have omitted the variable z throughout. Note 
that in going from the first line of (4.34) to the second line, we have used (4.12). 
Our next task is to interpret the numerator of the right hand side of (4.34), which 
is necessary because performing a direct asymptotic analysis of this quantity by 
using Propositions 4.1 and 4.2 proves to be difficult. That this is so is due to 
the degeneracy of the tridiagonal matrices £, and M. (Note that neither £ nor its 
factors C and M satisfy the genericity assumption in [MM].) For this purpose, we 
introduce for each 1 < j < ^' — 2 the shifted matrix £^^^ whose {k, I) entry is given 

by 

{£^^% = £k+j,i+j. (4.35) 
On the other hand, let '(/''■'' denote the solution of 

£^^u = zu, u-i =0,uo = 1. (4.36) 

Proposition 4.8. For each < j < p — 2, 

where 



<^j+i(2) V'j+i(^). 
_ forj odd (4-38) 

Proof. Prom the definition of £^^'^ and ^I-'^, it is clear that 

V't'^ {z) = ci {z)4>k+j [z) + C2 {z),Pk+j {z) (4.39) 
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for some 01(2;) and C2{z). By imposing the initial conditions in (4.36), we find that 

Therefore, on substituting into (4.39), we obtain 

Bj{z)i^^^\z) = -il^j.^{z)ct>k+j{z) + ^j.i{z)iPk+j{^). (4.41) 

Hence (4.37) follows from (4.34) if we replace j by j + 1 and let = p — 1 in 
(4.41). To complete the proof, it remains to calculate Bjj^-i{z). Here we make use 
of the quantity Wj{z) introduced in the proof of Proposition 4.3 which is related to 
Bj^i(z) by the relation Wj{z) = pjBj^i{z). From the proof of Proposition 4.3 (see 
relations (4.14) and (4.15)), we learn that W2j-i{z) is independent of the value of 
j, and the same holds true for W2j{z). Consequently, we have 

p,B^+,{z)=p,_2Bj-,{z). (4.42) 

Prom this, we find 



PjBj+i{z) 



P-iBo{z), for j odd, 
PqBi{z), for j even, 
Pp-i, for j odd, 
— ^^2^ , for j even. 



#(.) = + . . . + ^22±2^zi5^^, 

P2j ■ ■ ■ P2j + 2k-l P2j ■ ■ ■ P2j+2k-l ^ 

V'?^li(-) = ^^^±H^.'= 1^, 

^ P2j ■ ■ ■ P2j+2k P2j ■ ■ ■ P2j+2k Z''+^ 

(h) For < j < p/2 - 1 and < k < p/2 - 1, 



P2j+1 ■ ■ ■ P2j+2k P2j+1 • • • P2j+2k 



,/,[2j + l]/x_ Q'Sj fc+l «2j+2fc+l 1 

'y2k+i \Z) — Z ■ ■ ■ , . 

P2j+1 ■ ■ ■ P2j+2k+l P2j+1 ■ ■ ■ P2j+2k+l ^ 



In particular, 



(4.43) 



This completes the proof. □ 

In our next result, we will analyze il^^^^{z). 
Proposition 4.9. (a) For < j < p/2 - 1 and0<k<p/2- 1, 



(4.44) 



(4.45) 



[2,]^^ = _ a2,-2P2,-l «2,-1^2,-l 1 
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and 

+1] (^) = _ ^^2^ _J_ . (4.47) 

Proof. For the equation S^^^u = zu associated with the shifted matrix £"['1, the 
recurrence relations in (4.3) have to be replaced by 

P2k-lU2k-l-i — Oi2k-lU2k-i = z{a2kU2k-i + P2kU2k+l-i) , . 

_ _ (4.48j 

"2fc+l^t2fe+l-j + P2k+lU2k+2-i = z{p2kU2k-i — Q!2fe'f^2fc+l-i)- 

For i = 2j, we can therefore obtain (4.44) from (4.9) by shifting the indices. For 
i = 2j + 1, we have to solve (4.48) and an inductive argument leads to (4.45). □ 

Proposition 4.10. For j = 0, - ■ ■ ,p/2 — 1, we have 



^ ^ (4.49) 



p/2-J-l 



as z — >■ oo. Hence f2j{P) o,nd f2j+i{P) have poles of order p/2 — j — 1 at P+. 
Proof. From (4.46), (4,47), and Proposition 4.7, we have 

/2^(^)/2^(-)--^, /2^ + l(-)/2i + l(-) - (4-50) 

as z ^ oo. Therefore the assertion follows from (4.50) and (4.30). □ 

We next investigate the behaviour of f^iz) as z ^ 0. 
Proposition 4.11. For j = 0, • • • ,p/2 — 1, we have 



f2,{z) ^ -a2,-i I n ft I 



p-2 

n 

i=2j 



(4.51) 



\i=2j + l J 

as z ^ 0. Hence at Q_, f2j{P) has a zero of order p/2 — j, while f2j+i{P) has a 
zero of order p/2 — j — 1. 

Proof. By (4.25), Proposition 4.1 and 4.2, we have 

'"-^2,(^)~-|ap-ipa2i-i I lift) ^"^'"^ (4-52) 



i^p-i{z) 



32 



L.-C. LI, I. NENCIU 



as z — 0. On the other hand, 

On combining (4.31) and (4.32) and simphfy, we obtain the first relation in (4.50). 
The other relation in (4.51) follows in the same way. □ 

Proposition 4.12. For j = 0, - ■ ■ ,p/2 — 1, we have 

1 



^ ' (4.54) 

as z — ^ 0. Hence at f2j{P) has a pole of order p/2 — j while f2j+i{P) has a 
pole of order p/2 — j — 1. 

Proof. Prom (4.46), (4.47) and Proposition 4.7, we find 

fUz)f,-^iz) ~ /2++i(^)/2^+i(^) ~ ^ (4.55) 

as z — 7- 0. Therefore the assertion follows from (4.55) and (4.51). □ 

Combining Propositons 4.5, 4.7, 4.10, 4.11 and 4.12, we obtain the main result 
of the section. 

Theorem 4.13. For < j < p/2 - 1, 

(fe)>-i)-(|-i-i)p+ + (|-j-i)p- 

'p 
.2 



(/2,+i) > - (I - J - i) P+ + (I - J - i) P- 

where D = 'YlFk=i Pk- 

Corollary 4.14. For < j < p/2 — I, 

{f2j+p) >-D + {j + 1)P+ - (j + 1)P_ + jQ+ - jQ_, 

(/2,+1+p) >-D + {j + 1)P+ - {j + 1)P_ + (j + 1)Q+ - {j + 1)Q. 



(4.56) 



(4.57) 



Proof. This follows from (4.51), the relation /fc+p = h ^ fk for all k and (4.26). □ 
To close, we present the following result which is essential in Section 6 below. 
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Proposition 4.15. For each < j < p/2 — 1, the divisors 
and 

"l = ^+(f-.-i)n + (f-i)«.-(f-i-i)p--(f-i)Q- 

are general, i.e., 

dim L{U() = dim L{ui) = 1 (4.58) 
where for a divisor lA on C, 

L{U) = {meromorphic function cp | {(p) > —U }. 

Proof. We will adapt an argument of [MM] to our situation. For a divisor U on 
C, denote by Q,{U) the set of meromorphic 1-forms a; on C such that (w) > U. 
Take k and k' such that k + k'>g — 2=p — 3, then D + kP- + k'Q- has degree 
g + k + k' > 2g — 2. If (uj) > D + kP- +k'Q-, then co must be since a holomorphic 
1-form can have at most 2g — 2 zeros. Thus dimQ{D + kP- + k'Q-) = 0. By 
Riemann-Roch, it follows that dimL(D + kP- + k'Q-) = k + k' + 1 > g — 1. For 
concreteness, take = fc' = |, and we claim that 

L{D + (j - 1)P- + (j - 1)Q_) C L{D + jP- + (j - 1)Q_) 

P (4.59) 

CLiD+jP-+jQ-), 1<J<^- 

To establish this claim, wc just have to observe that by Corollary 4.14 above, we 
have /2(j_2)+i+p G L{D + (j - 1)P_ + (j - 1)Q_), but /2(,_2)+i+p ^ ^(Z? + 
jP- + {j - l)Q-). Similarly, /20-i)+p G L{D + jP- + (j - 1)Q_), but /2(i-i)+p ^ 
L(D + jP_ + jQ-). Thus it follows from dimL(D + f P_ + |Q_) = p + 1 and the 
claim in (4.58) that dim L(D) = 1. 

We next show that dimL{D + Q^ — Q-) = 1. Here we use the fact that allowing 
an extra pole increases dimL(ZY) by at most one. Thus we have 

1 < dimL(L> + (5+ - Q-) < dim L{D - Q-) + 1. (4.60) 

But L{D - Q-) C L{D), as /p_i = 1 G L{D) but /p_i ^ L{D - Q-). Hence we 
conclude from dimL(D) = 1 that dimL[D — Q-) = 0. Consequently, it follows from 
(4.54) that dimL(Z) + Q^ — Q-) = 1. Based on this, we can establish dimL{D + 
Q+ — Q- + P+ — P_) = 1 from the inequality 

1 < dimL(P> + Q+-Q-+P+- P_) < dimL(i:> + Q+ - Q_ - P_) + 1 (4.61) 

and the observation that L{D + Q^-Q-- P_) C L{D + Q+- Q-). (This follows 
because /p_2 G L{D + Q+ - Q-) but /p_2 ^ L{D + Q+ - Q- - P-).) Proceed 
inductively, we have the assertion. 

□ 
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5. Action-angle variables. 

As we saw in Section 2 above, the periodic Ablowitz-Ladik equation can be 
expressed in Lax pair form with Lax operator given by £{h). Therefore, the char- 
acteristic polynomial det(2;/ — £{h)) is invariant under the Hamiltonian flow and 
provides us with a collection of conserved quantities. Prom [S2], we have 

det(z/ - £{h)) = \ ]\pj\z^ [A(z) -{h + h'^)] 



(5.1) 



where the functions Ij as defined in the second line of (5.1) are such that 

= = 1, Ij=I-3^ i = 0,--- ,p/2-l. (5.2) 

Moreover, they are all polynomials in the a^'s, their conjugates, and P = 11^=0 Pji 
J =0,--- ,p-l. 

By using the fact that the collection ofpxp Floquet CMV matrices is a symplectic 
leaf of the Sklyanin bracket {•,-}j», we begin by reproving the involution theorem 
in [N] and [S2]. 

Theorem 5.1. The functions P, Iq, Relj, Imlj,j = I,-- - ,p/2 — 1 provide a col- 
lection of p conserved quantities in involution for the Ablowitz-Ladik equation. 

Proof. Write 

p 

det{zl - £{h)) = Er{£{h))zP-\ (5.3) 

r=0 

Then up to signs, the ^^^'s are the elementary symmetric functions. Prom (5.1) and 
(5.3), we find that 



f dh 
Ij{£) = j> _ Ep/^, .(£(/,))__, j = 0,... ,p/2 - 1, 

(5.4) 



'\h\ 

Hence the functions P,Io,ReIj,ImIj,j = l,--- ,p/2 — l are the pullbacks of central 
functions on to the 2p dimensional dressing orbit consisting of p x p Ploquet 
CMV matrices. Consequently, the assertion follows from the abstract involution 
theorem in [STS2] and Theorem 3.6. □ 

Remark 5.2. (a) Por readers who are not familiar with the abstract involution 
theorem in [STS2] , let us remark that its demonstration is a one line proof making 
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use of the fact that for a central function, its left gradient is the same as its right 
gradient. (See the expression for the Poisson bracket {-, -Ijh in (3.16) above for our 

case.) 

(b) Note that we are using the symbol £^ as a shorthand for the unitary loop £{■) in 
(5.4) above and we will henceforth continue to use the symbol with this meaning. 
Since wc will not be using the extended CMV matrices in what follows, this should 

not cause any confusion. 

Thus the number of commuting integrals as provided by the quantities in the 
above theorem is exactly equal to one half the dimension of the phase space. In the 
rest of the section, we will construct the variables (essentially) conjugate to these 
actions. As in Section 4, we denote by C the hyperelliptic Riemann surface of genus 
g = p — 1 corresponding to the affine curve I{h, z) = hdet{zl — £{h)) = 0. On C, 
we introduce the holomorphic 1-forms 

^k = -Q^dz, k = l,--- ,g = p-l. (5.5) 

dh 

We also introduce the meromorphic 1-form 

^p/2-l 
dh 

with poles at P± and Q±. Pick a fix point Pq on the finite part of C and put 
Dq = qPq. Then for £{h) satisfying the genericity assumptions {GA)i and (GA)2, 
we define 




where D = ^f.^^ Pk is the divisor of poles of /j, < j < — 1 in the finite part 
of C. Note that in the definition of v{£)-, the paths of integration going from the 
points of Dq to the points of D must avoid the points P±. These multi- valued 
variables are well defined because the points in Dq and D are in the finite part of 
C. On the other hand, the multi-valuedness can be resolved in the standard way 
and we will not try to get into the details here. (See, for example, [DLT] and [L2].) 
To compute the Poisson brackets between the conserved quantities in Theorem 5.1 
and the variables in (5.7), we will make use of a device in [DLT] (which has also 
proved to be successful in [L2]) which will allow us to simplify the calculation. In 
the following, we will deal with £{h) for h not necessarily on the unit circle. Note 
that in this general case, we have the relation 

£{h)£{h-^)* = £{h-^y£{h) = I, heC\{0} (5.8) 

which can be checked by using the fact that £{h) is unitary for h € 9D. For our 
purpose, we pick a fixed Hq £ (—1, 1) \ {0}, zq € dD such that (/iq, ^^o) is not on C 
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(5.9) 



and define 

Hho,zo{£) = Relogdet {zqI -£{ho)), 

Jh„,zo{£) = Imlogdet {zqI - S{ho)). 
As the reader will see in the calculation which follow, this choice of ho and Zq is 
critical. 

Lemma 5.3. (a) The Hamiltonian equation generated by Hh^^zoi^) is given by the 
equation 

£{h) = [£{h),B{h)] (5.10) 

where 

B{h) = {i{zoI - £{ho)r'£{ho)) _+iIm {i{zoI - £{ho))-'£iho))^ 

{i{zo - £{ho))-^£iho)) _] * + izoizol - f (V • ^^^^^ 

(b) The Hamiltonian equation generated by Jho,zo{£) is given by the equation 

£{h) = [£{h),C{h)] (5.12) 

where 

C{h) = {{zol - £{ho))-^£{ho)) _+iIm {{zol - £{ho))-''£{ho))o 

(5.13) 



{{zol - {ho))-'£{ho))_ - zo{zoI - £{ho')) 



1 

1 — hho 



Proof, (a) If we let (j){g) = Relogdet(zo-f — 5(^0)) for g G G^, then from the second 
equation in (3.60) and a direct calculation, we find that the Hamiltonian equation 
generated by Hh„,zo is given by (5.10), where 

B{h)=U~ ^\^^-£iho))-'£{ho)h 



h — h 







--{i{zoI - £{ho))-^£{ho))o + ilm (i(zo/ - £ (ho))-^ £ (ho)) ^ 
{i{zoI - £{ho))-'£{ho))_] * + (i{zoI - £{ho))-'£{ho 



h 



h — hi 





(5.14) 

Note that in going from the first line of (5.14) to the second line, we have used (3.12), 
together with the formula for Ilf, from [LI]. In the next step of the calculation, we 
will try to rewrite the last term in the above expression in the desired form, and it 
is here that the choice of ho and zo is important. To wit, by using (5.8), we have 

'i{zoI-£{ho))-^£{ho) ^ 



h — h, 



i{zo£{ho') - 1)-' 



h 



(5.15) 



h — ho 
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Hence (5.11) follows from (5.14) and (5.15). 

(b) The proof is similar to (a) and so we will skip the details. 

□ 

In our next two results, we will denote the Poisson bracket on the set oi p x p 
Floquet CMV matrices induced from {•, -Ijh simply by {•,•}. 

Proposition 5.4. For j, k = 1, ■ ■ ■ ,p/2 — 1, we have the following Poisson bracket 
relations 
(a) 

(b) 

(c) 
(d) 



lo, -Im (%^) } {S) = 1, {lo, -Im (i^±^^ } {£) = 0, 
)}(^) = 0; 



i>k—4>p-k 



Relj-Im (%^) } {8) = 0, {Uelj, -Im (^^±|i^) } (S) = <5,- ^/s-fe, 
Relj, -Re (^^iz|£^^ | (^) = Q; 

/m/,-, -/m (%^) } (^) = 0, |/m/,, -/m } (^) = 0, 

Iml„ -Re } {£) = 5,^p/2-k; 

P, -Im (%^) } {£) = 0, {P, -Im } (S) = 0, 

P,-i?e(^^^|2^)} (f) = 0. 

Proof. In order to compute }(£^), it suffices to evaluate it on an open 

dense set consisting of Floquet CMV matrices £{h) for which 

(a) the points Pj of the divisor D are distinct, 

(b) suppZ)n{f = o} = 0, 

(c) {QjVj=i n suppD = 0, where Qj = {hQ\zj{ho)) gC, j = !,■■■ ,P, 

(d) {Q,}f=in{f = O} = 0. 

So suppose £{h) satisfies (a)-(d) above. Then in the neighborhood of each Pj, 
we can take z to be the local coordinate and express h in terms of z. Thus in 
particular, hj = h{zj). Let £{t) be the Hamiltonian flow generated by IIho,zo 
let D{t) = Yl^Zi Pj{^) (where t is small) be the divisor of poles in the finite part of 
C of the corresponding eigenvector with last component normalized to 1, Pj{t) = 
(/i,(t),z,(t)). Then 



;=i dhy''J 



5 ^3) 



dt 



|t=o 



To compute the rate of change of Zj{t) at t = 0, consider an eigenvector f{z,t) 
for z in a neighborhood of Zj such that (ep_i, f{zj{t), t)) = for small values of t. 
Differentiate this relation with respect to t at t = 0, wc obtain 

dzj{t) _ (ep_i,5(/i,)/(z„0)) 



dt 



t=o (ep_,,|£(^^.,o)) 
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Therefore, on substituting this expression into (5.16), we find 

Now, on using the fact that the last column of a strictly lower triangular matrix is 
the zero vector, and {ep-i,R.e(^i{zoI — £{ho))~^£{ho))of{zj,0)) = 0, it follows from 
(5.11) that 

Bp-i, (^i{zoI - S{ho))~^£{ho) + izo{zoI - £{hQ^))~^ ■ ^ ^ f{zj,0) 

(5.18) 

Consequently, when we substitute (5.18) into (5.17), the result is 

p—i p—i 
{//^,,,„,</.fc}(^) = i^Resp^.Ffc + i^Resp^G'fc, (5.19) 
j=i j=i 

where 

Fk = {ep-i,{zoI - Siho))-^£{ho)v{z,0))Ck, 

Gk = {ep-i,zo{zoI - £ihQ^))~'^v{z,0)) ^ ^ ^ 

are meromorphic 1-forms on C. By a similar calculation, we also have 

p—i p—^ 
{Jho,zo,<t>k}{£) = ^Rcsp^Ffe - ^Rcsp^.Gfc. (5.21) 

Now from [S2], we know that for \zq\ = l,A(zo) € M, and so this implies Jho,zo 
is a constant independent of £(h). Consequently, { Jha.zo-, 4>k = and hence it 
follows from (5.21) that Resp^.F^ = Resp^G^. Consequently, 

p-i 

{Hh„,„<Pk}i£) = 2i^Resp^Gk. (5.22) 

Consider the meromorphic 1-form Gk- Since (|^)^ ~ ^1 as z ^ 0, it follows that 
~ ^z^~^ dz z ^ 0. On the other hand, we have 

^ -^2^/2^ ^1 as z ^ 0. (5.23) 



1 — h^h^ ho ' 1 — hoh 

Since the /^ 's are analytic at Q-, it follows that Gk has no poles at Q_. At the 
most singular component of ^(z, 0) is foiz) ~ trip-^ — \Z~'^^^. Hence it follows 

'-'p-l (,1 li = P'- ) 

from the asymptotics of foiz) and (5.23) that 

G^ ~ const, z''-^ as z ^ 0. (5.24) 
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So again Gk has no poles at (5+ . Thus Gk has poles only at D and at , j = 1, • • • , p 
as wc can similarly check that there are no poles at P±. Hence by the residue 
theorem and (5.22), 

p 

{Hu,,,,AkW) = -2i^ResQ,Gfe. (5.25) 



To simplify notation, let 6^- ' = ResQ^Gk and put b^'^^ = bj . We first calcu- 



late b'^!^\ We have 



1 



^ ho zo-Zj{ho) §ji{ho^,Zj{ho)) 

1 Zo 



zjjho)''-^ 

ho zo-Zj{ho) %{ho^,Zj{ho)) 



(5.26) 



Therefore, 



6^=) = J2 



Zj {ho 



ifc-i 



f^^ho Zo- Zj{ho) %{hQ ^,Zj{ho)) 



Zo_ 

ho 



lim 



dz ^ z, 



fc-i 





°^J\z\=R {zo - z) ■ I{ho ^ , z) 27ri /(/i^ ^,zo) 



(5.27) 



/io/(/io \ -2o) det{zoI - S{ho)) 
and so finally we conclude that 

{/f,,„,,„,</>fe}(^) = -2i6« 



-2z4 



det {zol - £{ho)) ■ 

But on the other hand, it follows from { Jho.zo^'Pk }{£) = and (5.1) that 



(5.28) 



{Hho,zo,<i>k }{£) 

--{Hho,zo +iJho,zoAk }{£) 

p/2 

{i,,<t^kW)zi-'''" ~ zi'\ho+h^'){PAkm I . 



det {zol - S{ho)) . ._ ,^ 

\3 — Pi ^ 

By equating (5.28) and (5.29), we conclude that 



(5.29) 



p/2 



-2iz^= ^ {I,-,</>fc}(f)z^+^/'-4/'(/io + /io'){^,'^fe}(^)> l<k<p-l. (5.30) 
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We now divide into three cases. 
Case 1: k>p/2 + l 

In this case, all brackets are zero except for { Ik-p/2-, (t>k }(^) = — 2z, i.e., 

}{E) = -2zJ,- fc_p/2, {Ij,4>k }{£) = 0, 1 < i < p/2 - 1, 
{P,</.fe}(f) = 0, {7o,0fc}(f) = O. 

Hence 

{Re/,-,(/)fc }(£:) = -i5j^k-p/2, {Im/j,^/)fe }(£^) = -Sj^k-p/2, 1 < j <p/2- 1, 
{P,<^fc}(5) =0, {/o,<^fc}(^:) =0. 

(5.32) 

Case 2: k<p/2-l 
In this case, we have 

{ Ij, ct>k }{£) = -2i5,- p/2-fc, { }{£) = 0, 1 < i < p/2 - 1, 
{P,(/.fc}(£:) = 0, {Jo,.^fe}(f) = 0, 

which implies 

{Rc/j,(/)fc }(f) = -i(5j,p/2-fc, {Im/j,^fc }(£:) = Sj,p/2-k, '^<j <pI2- 1, 
{P,<^fc}(^) = 0, {/o,<^fc}(^) = 0. 

Case 3: = p/2 

In this case, all brackets are zero except for 

-2i = {IoAp/2 W) - {ho + ho'){P,4>p/2 }{£)■ (5.35) 

Therefore, 

{ReIj,<Pp/2}{£)=0, {Im/„0p/2}(O=O, l<j<p/2-l, 

(5.3o) 

{P,4>p/2}{£) = 0, {loAp/2}{£) = -2i. 

□ 

Proposition 5.5. With v defined as in relation (5.7), we have the following Pois- 
son bracket relation: 

{p,Mg)}(^) = l. 

Proof. It suffices to compute the Poisson brackets { i^/io,zo , }{£) and { Jho,za , ^ }(^) 
on an open dense set consisting of Floquet CMV matrices £{}i) which satisfy con- 
ditions (a)-(d) in Proposition 5.4. Indeed, by following the same method of calcu- 
lation, we find 
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and 



where 



p—i p~i 
{ Jh,,,„i^ }{£) = ^ Resp.F - Resp.G = 0, (5.38) 



( -1 \ 



(5.39) 



Thus 



p-i 

{Hn„,,,v}{£) = 2i^Resp.G. (5.40) 
i=i 

Consider the meromorphic 1-form G. As z — >■ 0, it follows from the asymptotics of 
(see Section 4) and (f^)"^ that 

[rt-T^^dz. (5.41) 

Since the /j's are analytic at it follows from (5.23) and (5.41) that G has a 
simple pole at Q- Indeed, it follows from Propositon 4.11 that 



zMzoI-£{K'))-') ^ ^ 



ResQ.G = — li^^^. (5.42) 

In a similar way, it follows from (5.23), (5.41) and Proposition 4.12 that G also has 
a simple pole at (5+ and 

ResQ^G = ^ . (5.43) 

Now consider the two points P± at infinity and let u = ^ be the local coordinate. 
Then from (5.23), (5.41) and Proposition 4.10, we have 

ftjiz,0)-^^^ ^ —^^uUu as^^O, (5.44) 
1 - hoh^ ho [[■=2j Pi 

while 

f2j+iiz,0) ~ -—=^2 du asu^ 0. (5.45) 

1 - hoh^ ho Ui=2j+i Pi 

Prom (5.44) and (5.45), we conclude that G is analytic at P+. Similarly, by making 

use of Proposition 4.7 and (5.23), (5.41), we find that G has a simple pole at P_ 
with 

zo({zoI-£{h^'))-') 



Resp_ G 



p-l,p-2 Pp-2 



zo({zoI-£{K')) ') ^ ^ 

\ y p—1, p—1 



"^-^ (5.46) 
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Since G obviously has poles at the points of D and at the points Qj, j = 1, - ■ ■ ,p, 
it follows by the residue theorem, (5.40), (5.42), (5.43) and (5.46) that 



+ ((zo/-f(V'))"') , 



Pp-1 

hoap-i 

Pp-2' 



,p-2 ap-2- 

(5.47) 

To compute the second term on the right-hand side of (5.47), we introduce the 
following ad-hoc notation: if ^4 is a p x p matrix and < ji < ■ ■ ■ < ji < p — 1, 
0<:ki<---<ki<p— 1 are two sets of indices, with I > 1, then we denote 
hy A[ji, . . . ,ji;ki,... ,ki] the submatrix obtained from A by deleting rows ji, 
and columns ki, . . . ,ki. For simplicity of notation, let 

M = zol - S{ho^) . 

Then the terms appearing on the right-hand side of (5.47) can be identified as 

1 _ det(M[0,p-l]) , _ det(M[p-2,y-l]) 

K-1,0- dit(M) ^"""^ )p-hP-2- d^i(M) • 

To proceed, we expand both det{M[0,p — 1]) and det(M[p — 2,p — 1]) along their 
0th columns, which makes all the minors appearing the calculation /iQ-independent, 
and so allows us to separate the /iQ-dependent terms from the /iQ-independent ones. 
This leads to 



where the /io-dependent term is 

(o.4yj 



+ ho ■ pp-2Pp-idet{M\p - 2,p - 1; 0,p - 1])^ 



and 1712 is /io-independent. While the minors appearing in m2 have a structure 
which cannot be easily simplified, the minors appearing in mi can be computed 
explicitely. 

Indeed, for each < j < | — 1, consider the 2x2 blocks 

^, ^ f -P2j-ia2j zo + a2j-ia2j\ A=( ~P'^i^W -p2jp2j+i 

^ \-P2j-ip2j a2j-iP2j J ^ \ZQ + a2ja2j+i a2jP2j+i 

Note Aj and Aj are, respectively, the left and right "halves" of the 2x4 blocks 
appearing in the extended matrix zqI — £. In particular, direct investigation shows 
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that M[0, 1; 0,p — 1] is a block bi-diagonal matrix, having the blocks Ai, ... , Ar_i 

on the diagonal and Ai, . . . , Az_2 on the upper diagonal, while M\p—2,p—l; 0,p— 1] 

has Ao,. . . ,Ap_2 on the diagonal and Ai,. . . , Ap_2 on the lower diagonal. This 
implies that 

f-l f-l ^_ p-2 

det(M[0, l;0,p- 1]) = J] det{Aj) = ]J zoP2j-iP2j = 4~' H Pi ' 

and 

det(M[p - 2,p - 1; 0,p - 1]) = [] det(i,) = [] ^op2iP2,+i = H ■ 

j=0 j=0 j=0 

Plugging these two expressions in (5.49) leads to a very simple expression: 

£-1 

mi = ^ ^-TTT • (^0 + 7-) • (5.50) 

det{zoI - S{ho^)) y hoJ ^ ^ 

On the other hand, proceeding exactly as in the proof of Proposition 5.4, we find 
that 

y ResQ . G = ° ^ . (5.51) 

^tt ' det{zoI - Siho)) ^ ' 

By combining (5.48), (5.50), and (5.51) into (5.47), and using (5.38), we conclude 
that 

{1^.0.0 + ^Jn...,-m = det(.o7-g(^o)) t " ~ ' 

By the analogue for u of relation (5.29), together with the fact that m2 is ^q- 
independent, we conclude that 

{P,u}{S)=Ai, 

which leads directly to our claim. 

□ 

Remark 5.6. (a) In the Proposition above, we were unable to obtain the Poisson 
brackets between the Tj 's and v. Nevertheless, the functional independence of the 
integrals follows by combining Proposition 5.4 and Proposition 5.5. 
(b) Proposition 5.4 shows that the Hamiltonian flow generated by P does not give 
rise to nontrivial motion on the Jacobi variety of C, this is the reason for the in- 
troduction of the algebro-geometric variable v. (See also Proposition 6.2 below in 
this connection.) Note that in particular, this means that the periodic defocusing 
Ablowitz-Ladik equation is not an algebraically completely integrable system in the 
sense of Adler and van Moerbeke [AMV] . For other examples of integrable systems 
involving spectral curves which are not algebraically integrable, we refer the reader 
to [DLT] and [L2]. 
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6. Solving the equations via factorization problems. 

In this section, we will solve the Hamiltonian equations of motion generated by 
the commuting integrals in Section 5 via factorization problems. 

We begin by writing down the Hamiltonian equations of motion by using Propo- 
sition 3.9, Theorem 3.6 and (5.4). To do so, for a map F : gl{p,C) — > C, we let 
VF(M) = {dF / dniij) denote its gradient. In what follows, we will regard Re Ij, 
Im Ij and P as functions of g'^, g°, where £{h) = g^g°{h), and £{h) = g°{h)g^. We 
will use the notation of Proposition 3.9 (b). As an example, if i7 = Im Ip/2-j, then 
the corresponding central function tp is given by ip{X) = Im Ej{X{h))-^f^. 

Proposition 6.1. The Hamiltonian equations of motion generated by H = ImIp/2-j, 
Re Ip/2-j J 0,'f'd P O'f^ given by 



g^ = {U~,jD^{S{hmg^ - g%U~,jD^{S{hm, 

g°{h) = {U,{D^{S{hmg"{h)-g"{h){U-,(D^{Sm), 



(6.1) 



where 



' S{h)V^Ej{£{h)), for H = ImIp/2-j 

Dip{£{h)) = I iS{h)V'^Ej{S{h)), for H = ReIp/2-j (6.2) 

^ -ih£{h)V'^Ep/2{£{h)), for H = P, 

and similarly for Dip{£{h)). 

In (6.1), the parameter h is on the unit circle, however, we will remove this 
restriction later on. (See (6.29) below.) Before solving these equations, let us spell 
out V^Ei (x) more explicitly. To do so, observe that V'^Ej{x) (see (5.3) above) 
obey the recursion relations 

V'^Ej+i{x) = xV'^Ej{x) - Ej{x)I, 0<j<p, (6.3) 

where by convention V^£^p+i(x) = 0. Since V'^Eo{x) = 0, by solving the recursion 
relations backwards, we obtain 

J-i 

V''Ej{x) = -J2Ej-i-i{x)x\ j = !,■■■, p. (6.4) 

i=0 

The equations of motion generated by P are the simplest to solve. Although 
we could easily write down the solutions of these equations without recourse to 
Proposition 6.1 above, however, we will do it by the proposition in order to achieve 
uniformity in our treatment. 

Proposition 6.2. The Hamiltonian equations of motion generated by P simplify 
to 

r = <7'Ai-A25' 
g''{h)=g%h)A2-A,g%h) 
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where 

Ai = iRe (^Ai(^^^^^ = diag{0,iP,0,iP,--- ,0,iP), (6.6) 

and 

k2 = iRe (a_i{^^^ = diag{iP,i),iP,0,--- ,iP,Q). (6.7) 
The solutions of (6.5) are therefore given by 

g'^it) = e-*^2(°)/(0)e*^^(°),5°(/i,i) = e-*^^(°)5°(/i,0)e*^^(°). (6.8) 

Hence 

aj{t) = a,(0)e^*^(°), pj{t) = p,(0), j = 0, • • • ,p - 1. (6.9) 
Proof. According to (6.2) and (6.3), 

Dip{S{h)) =ihY^ Ej{£{h))£{h)^-^ . (6.10) 

3=0 

For j = 1, • • • , I - 1, it follows from (2.11) that 

(ihE,{£{h))£{hY^-^) = ihEj{£{h))£{hy^-^. (6.11) 
Similarly, by using the fact that A-i (|) is lower triangular, we find 

(jhEoi£{h))£{h)'^^ = ih£{h)'^ - ilm (iA_^ ^) ^ (^"^^^ 
as Eq = 1. Therefore, on using (6.11) and (6.12), we obtain 

Il~^^Dip{£{h))=D^{£{h))-iKe (6,13) 
By a similar calculation, we find 

Jl-,^D^{£{h)) = Dipim) - iRe (^1 (1))^ . (6.14) 

Hence, on substituting (6.13) and (6.14) into (6.1), and using the obvious facts that 
£{h)2-3g^ = g^£{h)2-j ^ g°{h)£{h)2-i = £[h)i~^ g°{h), we obtain the equations 
in (6.5). The formulas in (6.6) and (6.7) for Ai and A2 in terms of P then follow 
from Lemma 2.2 and (2.19). As P is a conserved quantity, it is easy to verify that 
the expressions in (6.8) give solutions to the equations in (6.5). Finally the solution 
formulas in (6.9) are obtained from (6.8) by multiplying out. □ 

To solve the Hamiltonian equations generated hy H = Im/p/2-j) R,e/p/2-j) we 
will make use of Proposition 3.9 (c), which means we have to solve explicitly for 
each t > the following factorization problems 

gtD<p(£(M)) = ki{h,t)h{h,t)-\ e'^'^^^^^'^^^ = k2{h,t)b2{h,t)-^ (6.15) 
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for ki{-,t) G Kyj, bi{-,t) G -B^,i = 1,2. However, from the definition of K^, it is 
easy to show that it is enough to solve 

^-tD<pisih,o))^-tD^(eih,o)r ^ bi{h,t)h{h,ty, 
^-tD<pieih,o))^-tD^ii(h,o)r ^ i^^^h, t)b2(h, ty ^^'^^^ 

for bi{-,t),b2{-, h) G 5^. Note that for i = 1,2, bi{-,t) (resp. can be ex- 

tended analytically in the interior (resp. exterior) of the unit circle \h\ = 1. So 
(6.16) is a Riemann-Hilbert problem. In order to solve this problem explicitly, we 
will first transform the product on the left hand side of (6.16) into a form which 
makes the problem more tractable. To this end, note that it follows from (6.4) that 

i-i 

V''Ej{£{h,0)) = -Y,Ej-i-i{S{h,0))Sih,Oy, j = l,---,p. (6.17) 

As £{h,0) commutes with £{h,Oy for h G dB, it follows from (6.17) and (6.2) that 
Dip{£{h, 0)) commutes with Dip{£{h, 0))* for h G 5D. In a similar way, we see that 
Dip{£{h,0)) commutes with Dip{£{h,0))* . Consequently, we can rewrite (6.16) as 

g-t(i?^(£(M))+/5v'(f(M))*) = bAh,t)bAh,ty, 

(6.18) 

^-t(D^i£ih,o))+D^isih,o)n = b^^h^ t)b2{h, ty, he do. 

Now if we compute D(p{£{h,0)y (resp. D(p{£{h,0))) more carefully, we find 

Dip{£{h, 0)y = ±Dip{£{h, 0)-^), 

Dip{£{h,(})y =±Dip{£{h,0)-^), hedB, 

where we pick the + sign for H = lmIp/2-j and the — sign for H = KeIp/2-j- 
Hence (6.18) becomes 

^-t{D^{S{h,0))±D^{S{h,0)-^)) ^ bi(h,t)bAh,ty, 

(6.20) 

g-t(D^(£(M))±o^(f(M)-^)) = e 

where the choice of sign is described in the previous sentence. 

As the explicit solution of the factorization problem in (6.20) will involve con- 
structing bi{h,t),i = 1,2, for values of h not on unit circle, it is necessary to 
introduce some Lie algebras and projection operators which complements those in 
Section 3. For this purpose, let A be the ring of Laurent polynomials in the vari- 
able h and let glp{A) be the Lie algebra of p x p matrix functions with entries 
in A equipped with the pointwise Lie bracket. We will consider the following Lie 
subalgebras of glp{A): 



^Xjh^ G glpiC[h]) \Xoeb 

3=0 

l=[X e glp{A) I X{h)+X{h-^y = 0} . 



(6.21) 
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Then analogous to (3.8), we have the sphtting 

glp{A) = teb. (6.22) 

For X e glp{A), define 

(n+x)(/i) = Y.x,h\ (n-x)(/i) = Y^x.h^, (6.23) 

j>0 j<0 

then the projection operator onto b associated with the sphtting in (6.22) is given 

by 

{U~^X){h) = iU+X){h) + HfaXo + {iU_X){h-^)r . (6.24) 

Theorem 6.3. For H = ReIp/2-j,ImIp/2-j , there exist unique holomorphic matrix- 
valued functions 

bi{-,t) : CP^ \{oo} — > GL{p,R), 
b2{;t) :CP^\{oo} — >GL{p,R) 



(6.25) 



which are smooth in t, solve the factorization problems 

^-t(D<p(sih,o))±D^ieih,o)-^)) ^ b^(^h,t)bi{h-\ty, 

g-t(D^(£(h,0))±D^(£>,0)-)) (6.26) 

heCF^\ {0, oo} 

(where the + sign corresponds to H = ImIp/2-j and the — sign corresponds to 
H = ReIp/2-j) and satisfy 

bi{0,t) e B,bi{h,t)~^bi{h,t) e ImU^ i = l,2. (6.27) 

Moreover, for h G CP^ \ {0, oo}, the formulas 

g^t) = br{h, t)-'g%0) b2{K t) = b^{h-\tr g<'mb2{h-\t)*)-\ 



g%h, t) = b2{h, t)-^g%h, 0) 6i {h, t) = b2{h-\t)*g%h, 0)(6i (h-^t)*)-^ 

give solutions of the equations 

r = {U~^D^{S{-mh)g^ - g%U~^D^{£{.mh), 

g%h) = iU~,D^{£{.mh)9°{h)-g%h)iU~,D^{Si-))){h). 



(6.28) 



(6.29) 



Finally, for generic initial data ^^(0) and g"{h,0), bi{h,t) and b2{h,t) can be con- 
structed by means of theta functions associated with the Riemann surface of the 
spectral curve C = {{h,z) \ det{zl — g'^{0)g°{h,0)) = 0} for values of t > for 
which Oij{t) / 0, J = 0, • • • ,p — 1. 
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Proof. We will prove the result for H = Im/p/2-j- The argument for the other case 
is similar. We start with uniqueness of the factors bi{h,t),i = 1,2. To prove this, 
suppose b'^{h,t),i = 1,2 is a second pair of solutions of the factorization problem. 
Then from hi{h.t)h{h-^ ,t)* = Wi{h,t)b1(h-^ ,t)* , we have 

gi{h,t) = b'l{h,t)-%{h,t) = b'l{h-\ty{bi{h-\ty)-\ /i G \ {0,oo}. (6.30) 

Clearly the function gi{-,t) defined in (6.30) above can be extended to an analytic 
function everywhere, hence by Liouville's theorem, gi{h,t) = Ci[t). To determine 
Ci{t), note that gi(0,t) = 6"(0, t)~^6i(0, t) G b. On the other hand, gi{oo,t) = 
{gi{0,t)~^)* is upper triangular with positive diagonal entries on the diagonal. 
Hence Ci{t) = /. 

To establish the existence of bi(h,t),i = 1,2, note that S{h,t) (resp. £{h,t)) 
exists for h G CP\{0, oo} since it exists for h G dB). Hence we can obtain bi{h, t),i = 
1, 2 as solutions of the equations 

k{h,t) = -b^{h,t){Ut,D^{S{;t))){h), (6.31) 



and 

b2{h,t) = -b2{h,t){UbDip{£{h,t))){h). (6.32) 

Clearly, the analyticity properties and (6.27) are satisfied by definition. We next 
consider the product bi{h,t)bi(h~^ ,t)* . By differentiating, we have 

= - h{h,t) {(IliD<p(£{;t))){h) + ((IliDy,{e(;t))){h-^)r) h{h-\tr. 

Now, from the definition of H^, it is straightforward to check that 

{U~,X){h) + {{U~,X){h-')r = X{h) + {X{h-')r. (6.34) 
Consequently, we obtain 

= Dip{£{h,t)) + {Dip{£{h-^,t))y (6.35) 
= Dip{£{h,t)) + Dip{£{h,t)-^). 

Substitution of (6.35) into (6.33) therefore gives the relation 

±(b,ih,t)b^ih-\tr) 

= -bi{h,t){Dip{£{h,t))+D^{£{h,t)-'))bi{h-\t)* (6-36) 
= -bi{h,t)bi{h-\ty {Dip{Tih,t)) + D(p{J^ih,t)-^)) 
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where we have used the fact that ip is a central function, and where F{h,t) = 
{bi{h-\t)*)-'^S{h,t)bi(h-^ ,t)* . Now, by direct differentiation, using (6.31), the 
equation 

Sih,t) = iU~^D^{£{-,tmh)£{h,t)-SihM'^~,Dip{£{-,tmh), (6.37) 
and (6.35), we find that 

j^T{h, t) ={hih-\ty)-\Dv{£ih, t)) + Dip{£ih, t)-^))S{h, t)h {h-\ty 

- {hCh-\ty)-^£{h,t){Dip{£{h,t)) + Dcp{£{h,t)-^))bi(h-\t)* 
=0. 

(6.38) 

Therefore, F{h,t) = £{h,0) and so (6.36) becomes 

= -bi{h,t)bi{h-\t)* {Dip{£{h,0)) + Dip{£{h,0-'))) . 

This shows bi{h, t) satisfies the first relation in (6.26). In a similar fashion, we can 
show that b2{h,t) satisfies the second relation in (6.26). 

Finally, we will show that g^{t) and g°{h,t) as defined in (6.28) satisfy (6.29). 
First, note that by using the relation £{h,0) = g^{0)~^£{h,0)g^{0) and (6.26), we 
have 

biih,t)-'g%0)b2{h,t) = biih-\tyg%Q){b2{h-\tr)-\ 

b2{h,t)-'g"{h,o)bi{h,t) = b2{h-\tyg"{h,oybi{h-\ty)-\ 

Differentiate ^^(t) = bi{h,t)-^g''{o)b2{h,t) and g°{h,t) = b2{h,t)-^g°{h,0)bi{h,t) 
with respect to t, we find 

g%t) = -b^{h,t)-%{h,t) g-{t) + g''^{t)b2{h,t)-%{h,t), 
g%h,t) = -b2{h,t)-%{h,t)g"{h,t)+g%h,t)bi{h,t)-%{h,t). 

On the other hand, by differentiating the first relation in (6.26) with respect to t, 
and multiply the resulting expression on the left by bi{h,t)~^ and on the right by 
(6l(/l-^^)*)-^ we obtain 

-b,Ch-\ty{D^{£{h,o)) + D^{£{h,o)-^)){b,Ch-\tyr' 
=bi{h,t)-%{h,t) + h{h-\ty{bi{h-\ty)-\ 

But from the fact that ip is a central function and (6.40), we see that 

b,{h-\ty{Dip{£{h,0)) + Dpi£{h,0)-')){b,{h-\ty)-' 

1 (6.43) 
= D(p{£{h,t))+Dip{£{h,t)-^) 
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where S{h,t) = g''{t)g"{h,t). Thus (6.42) becomes 

-{Dip{£{h,t)) + Dip{S{h,t)-^)) = h{h,t)-^h{h,t) + {h{h~\t)-^h{h-\t))*. 

(6.44) 

Now let 

xih,t) = hih,t)-%ih,t) - ihih-\t)-%ih-\t)r, 

Y{h,t) = Dip{£{h,t)) - Dip{£(h-\t)-^). 

Clearly, X{-,t) G k, thus U-^X{-,t) = 0. On the other hand, as Dip{£{h,t)-'^) = 
D(p{£{h~^ ,t))* , we also have Y{-,t) G k and Il-^Y{-,t) = 0. Consequently, when we 
apply to both sides of (6.44), we obtain 

h{h,t)-%{h,t) = -{U^Dip{£{;t))){h). (6.46) 

Similarly, from the second relation in (6.26), we can show that 

b2{h,t)-%{h,t) = -{U^Dip{£{;tmh). (6.47) 

Finally, substituting (6.46) and (6.47) into (6.41) gives (6.29). This completes the 
proof of the theorem modulo the assertion on the construction oibi{h,t) and 62 (/i, t) 
via Riemann theta functions. □ 

We now turn to the construction of bi{h,t) and b2{h,t) via theta functions. 
Again, we will give details for H = lmIp/2-j, leaving the other case to the interested 
reader. The following proposition shows we can construct b2{h,t) from bi{h,t) and 
the solution of a finite dimensional factorization problem. 

Proposition 6.4. Let l{t) £ B, u(t) G K be the solution of the factorization 
problem 

biiQ,t)-^g%0)=u{t)l{t)-\ (6.48) 

Then 

b2{h,t) = g%Orbi{h,t)uit). (6.49) 

Proof. Since £{h,0) = g^{0)~^£{h,0)g^{0), the factorization problem for b2{h,t) in 
(6.26) can be rewritten as 

^-t{D^{£{h,o))+D^{e{h,o)-')) ^ g-(^o)b2{h,t)b2{h-\t)*{g%0))-\ (6.50) 
Therefore, when we compare this with the first relation in (6.26), we obtain 

b2{h,t)b2{h-\ty = g%oybi{h,t){g^{orbi{h-\t)y. (6.51) 

Now let l{t) G B, u{t) E K he the solution of the factorization problem in (6.48). 
Then 

/(0)*6i(/i, t) = m + Oih))u{ty (6.52) 
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and 

g%0)*h{h-\t) = {l{t) + 0{h-^))u{t)* . (6.53) 
Substitute (6.52) and (6.53) into (6.51), we obtain 

b2ih, t)b2{h-\ t) = m + Oih)){l{t)* + 0(/i-i)). (6.54) 

Therefore, from the uniqueness of solution of the factorization problem (Theorem 
6.3), we conclude that 

b2{h,t) = l{t) + 0{h) = g%oybi{h,t)u{t). 

□ 

To construct bi{h,t), we invoke the formula in (6.2) and (6.17), according to 
which we have 

{D<p{S{h{P), 0)) + D^{S{h{P), 0)-^))v{P) = HiPMP) (6-55) 

where fij is meromorphic on the hyperelliptic Riemann surface C. Prom the first 
relation in (6.26), 

e-'^'^^^^bi{h{P),t)-^v{P) 
=5i (hiP) , t) -1 e-*(«^(^(MP),o))+Dv(£(/.(P),o)- ))^p^ (g_5g^ 



=bi{h{P) ,trviP) 
for h{P) e CP^ \ {0,oo}. Since S{h,t) = bi{h,t)-^£{h,0)bi{h,t), we have 

£{h{P),t)biih{P),t)-'v{P) = z{P)bi{h{P),t)-^v{P). (6.57) 
On the other hand, as 



s{h{P),t) ={£{h{P) ,try 

=biih(P,tr£ih{P),0)ibr{h(P)~\ty)-\ 

we also have 



(6.58) 



S{h{P),t)b^{h{P) ,t)*v{P) = z{P)bi{h{P) ,tyv{P). (6.59) 
Thus if we let 



v'^{P) = b^ih{P),t)-'v{P), 
vl{P) = b^ih(P)~\trv{P), 
then (6.56), (6.57) and (6.59) give 

e-tM,(P)^^(p) = /j(p) g cpi \ {o,oo} 

£{h{P),t)vi{P) = z{P)vi{P). 



(6.60) 



(6.61) 
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In this way, we axe led to scalar factorization problems. Note that because 6i(0,t) 
is lower triangular, it follows from (4.50) and (6.60) that 

(«)2,) > -D{0) + (I - J- - l) P_ + (I - Q_ 

(«)2,+i) > -D{0) + (I - J- - l) P_ + (I - J - l) Q_ 

on C \ {h{P) = oo}. Similarly, because bi{0,t)* is upper triangular, we find that 

{{vl)2,+i) > -D{0) - (I - J - l) P+ - (I - J - l) 
on C\{h{P) = 0}. 

We will first solve the following scalar factorization problem (cf. [RSTS], [DL]) 
g-tM,(p)^^(P) =uj'_{P), PeC\ {h{P) = 0,oo}, 
(w+) > -^(0), on C \ {h{P) = oo}, (6.64) 
(wi) > -D{0), on C \ {h{P) = 0}. 

To do so, we fix a canonical homology basis {aj, bk}i<j^k<g of the Riemann surface 
associated with the spectral curve, and let {wj}i<j<c, be a cohomology basis dual 
to {aj,bk}, i.e., 

/ Ui = 5ij, / uji = Qij, (6.65) 

where (^iy ) is the Riemann matrix. Let 

0{z) = 9{z,i},) = ^ exp{27ri(m, z) + '!ri{m,ftm)} (6.66) 

be the Riemann theta function associated with the matrix $7. Let uj = {oji, ■ ■ ■ ,ujg). 
Choose a nonsingular e G in the theta divisor, i.e. 0{e) = 0, the prime form 
Ee{x,y) = ^ (e + /Jw) ^0 with the additional property that Ee{P±,P) are not 
identically zero in P. (Sec Lemma 3.3 of [M].) Let Pq be a fixed point on the finite 
part of the Riemann surface, then by Corollary 3.6 of [M], there exists an effective 
divisor Dg-i of degree g — I such that e = A — /(^^j^jp^ where A is the vector of 
Riemann constants. Also note that 




r.(g-l)Po + P \ 

e+ / uj\=9[A+l uj\ =0 ^ P e D{0). (6.67) 

>(0) / 



Now let d(f>+ be the unique meromorphic differential of the second kind with van- 
ishing a-periods with poles only at {h{P) = oo} such that {d(f)+ —dfij){P) is regular 
in C \ {h{P) = 0}. Set 

c^P) = exp (tMP)) ^ . ^[j, +p X ^ (6-68) 
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where 4>+{P) = Jp^ d(f)-^- and is the vector of 6-periods of i.e., 

By using (6.69), the basic property 

e{z + y + O7) = 9{z) exp (^27ri(-(7, z) - ^(7, Q7)) j , 

and (6.67), it is straightforward to check that ui\_ is single- valued and meromorphic 
in C \ {h{P) = 00} with (oj^j^) > —-0(0) there. Moreover, e~*^j(^)cj^_(P) is mero- 
morphic in C\{h{P) = 0} with (6"*^^^^) > -D{0) because d(t)+{P) - dfij{P) is 
regular in C\{h{P) = 0}. Thus and uj*_ = e~*^^ solves the scalar factorization 
problem (6.64). 

Next, we compare (6.64) and (6.61), this gives 

coXiPy^v+iP) ='^-(P)~^v'_{P), C\{h{P) =0,00}. (6.70) 

TliGrGfor6 

?y*fpi ^ I P^c\{h{p) = oo} 

^ ' \ Lo'_{P)-^v'_{P), P G C \ {/i(P) = 0} ^' ' 

is meromorphic on C. Moreover, it follows from (6.62), (6.63) and the expressions 
for oj^. that 

(5*2,) > -m - (I - J - 1) p+ + (I - i - 1) p. 



(55,+i) > -m - (I - J - 1) ^+ + (f - J - 1) ^- 
-(f-i-i)Q..(i-i-0«- 

where -D(i) is the divisor of zeros of the function 



(6.72) 



u + t^+\. (6.73) 
(0) / 



Since 8{h{P),t)v*^{P) = z{P)v*^{P), it follows from the definition of if*(P) that 

S{h{P),tyS'{P) = z{P)^{P). (6.74) 

But on the other hand, 

S{h{P), t)v{t, P) = z{P)v{t, P), (6.75) 
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(6.76) 



where the last component of v{t, P) is equal to 1 and 

(»2j(«, ■)) > -0{t) - (I - i - l) P+ + (I - J - l) P- 

- (I <3-, 

(«,•)) > -m - (I - J - 1) ^+ + (I - J - 1) ^- 

~(f-.^-i)<3..(f-.-i)Q-. 

Clearly, v*{P) = v*p_i{P)v{t, P). Let DQ{t) be the divisor of zeros of v^-i{P) so 
that _ _ 

{%_,) = Doit) -Dit). (6.77) 

Then it follows from the relation connecting v*iP) and vit,P) above, (6.72) and 
(6.77) that 

iv2jit,-)) =i^2j)-i%-l) 

> - Doit) -(l-j-l)p^ + (l-j- l) P_ 



Similarly, 

(^^2,+i(t,-)) > -Doit) -(^-j-l)P+ + (^-j- l] p. 

(6.79) 



(f-,-i)«..(i-.-i)«- 



Clearly, we must have Doit) > -D(t). Since deg Dit) = deg -Do(i) = 9, we must 
have Doit) = Dit). But -D(t) is a general divisor by Proposition 4.12, so from 
iv'p-i) = Doit) — Dit) = Dit) — Dit), we must have S^.i = constant and this 
implies -D(i) = Dit). Thus we can solve for v^iP) and hence f±(-P) up to multiples 
because of Proposition 4.12. Indeed, by making use of the prime form, we can write 
down the explicit expression 

'E,iP_,P)\^-'-^ fE,iQ_,P)\^-' 



^2,iP) -cAt) y^J^)) • \-EjQ^)^ 

rD,-.+P+i%-3-i)P- + i%-o)Q- \ (6.80) 

^l^^ + JDW + (f-,-l)P+ + (f-j)Q+ 



Dg-i+P . 

(t) 



where C2j{t) has to be determined. Similarly, we have 

-t (,^(EAEi:iiyV~'~' f E,iQ.,P) \^-^-' 

V2,^.ip) =c2,^.it) [e:^^) ■ [e:(q^^) 

f) ^ .i?.-i+P+(f -j-i)P-+(f -j-i)Q- \ (6.81) 

^ \^ + ■/i3(t)+(§-,-l)p^ + (|-,-l)Q^ 
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where C2j+i {t) is also as yet undetermined. 

We are now ready to construct bi{h, t) and in the process, we will also determine 
Cj{t). For given h G CP^ which is not a branch point of the coordinate function 
h{P), there exist p points Po{h), - ■ ■ , Pp-i{h) of the Riemann surface C lying over 
h. Therefore we can define the matrices 

v±{h,t) = {viiPoih)) ■ ■ ■ vi{p,_,m 
v\h) = {v{Po{h))...v{p,.rm 

where v{P) can be obtained from the formula for v^{P) by setting t = (since 
a;jr°(P) = 1 and bi{h,t = 0) = /) and so can be computed in terms of theta 
functions. Of course, v^{P) are also in terms of theta functions. With these 
matrices, it follows that 

bi{h,t) = V\h)V+{h,ty^ . (6.83) 

Of course, v±(-P) = a;^(P)5*(P) are determined only up to the Cj(t)'s. Write 

viiP) = cit)vi{t,P) (6.84) 

where Vj.{t, P) are known and c{t) = diag (co(t), ■ ■ ■ , Cp_i(t)) is to be determined. 
Then 

V±{h,t) = c{t)Vl{h,t) (6.85) 
where Vl{h, t) = (v^(t, Po{h)) ■ ■ ■ vi{t, Pp-i{h))). With these definitions, 

bi{h,t) = V\h)Vfih,t)-^c{t)-\ (6.86) 
As bi{h,0) = I, the above relation determines c(0) via the formula 

c{0) = V\h)Vfih,0)-\ (6.87) 
To determine c{t) for t > 0, we bring in 

h{h,t) =e*-°*'(^(^'°»6i(/i,t) 

which is required to be unitary for h G dH). Therefore, if we equate the expression 
for /ci(l,t) from (6.88) with the corresponding one for {ki{l,t)*)~^ , we obtain 
|c(t)|2 = c{t)*c{t). Explicitly, 

|c(t)|2 = (y^(i)yf(i,t)-i)*e*(^'^(^(i'0))+^^(^(i'0)"'))y^yf(i,t)-i. (6.89) 

Write c{t) = |c(t)|e*''(*) , where e*''^*) = diag (e*''"^*), • • • , e*''p-i(*)). It remains to 
determine e"'^*-* . However, this is fixed by the condition that 

6i(0,t) = y^(0)Ff(0,t)-i|c(t)rie-^''W G B (6.90) 

as the diagonal entries of the elements in B are positive. 
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